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SUMMARY 


The investigation into the appearance of intruder states from the neg- 
ative continuum when some of the two-electron integrals involving the small 
component of the wave function were omitted has been completed and published 
[1] . The work shows that provided all integrals involving core contracted 
functions in an atomic general contraction axe included, or that the core 
functions are radially localized, meaningful results axe obtained and in- 
truder states do not appear. Reprints of the paper are attached. 

In the area of program development, the Dirac-Hartree-Fock (DHF) pro- 
gram for closed-shell polyatomic molecules has been extended to permit 
Kramers-restrict ed open-shell DHF calculations with one electron in an open 
shell or one hole in a closed shell, or state-averaged DHF calculations over 
several particle or hole doublet states. Code for two open shells is still 
being tested, as is code for limited MCSCF, for situations where the ground 
state at the DHF level cannot be described by a single determinant, but is a 
linear combination of a few determinants related by double excitations from 
one Kramers pair into another. One application of the open-shell code was to 
the KO molecule [2] . Preprints of this paper are attached. 

Another major area of program development which is under development is 
the transformation of integrals from the scalar basis in which they are gener- 
ated to the 2-spinor basis employed in parts of the DHF program, and thence to 
supermatrix form. The reason for these developments is that the DHF program, 
while written efficiently for the circumstances under which it was expected 
to be used, is not the most efficient implementation possible. The code was 



originally developed under the assumption that disk space would not be avail- 
able to sort and transform the integrals, and therefore the unordered inte- 
grals would have to be used in the construction of the Fock matrix. With the 
experience gained over the past year, particularly concerning the omission 
of small component integrals, and with the increase in availability of disk 
space, it is now possible to consider transforming the integrals. The use of 
ordered integrals, either in the scalar basis or in the 2-spinor basis, would 
considerably speed up the construction of the Fock matrix, and even more so if 
supermatrices were constructed. Furthermore, in order to proceed beyond the 
SCF level and include electron correlation it is necessary to transform the 
integrals to the molecular 4-spinor basis. Therefore, a considerable amount 
of effort has been spent on analyzing the integral ordering and transformation 
for the DHF problem. Much of this work was used in preparation for the NATO Ad- 
vanced Summer Institute in Vancouver, BC, in August 1992. 

The work of assessing the reliability of the relativistic effective core 
potentials (RECPs) available in the literature has been continued with calcu- 
lations on the group IV monoxides. The perturbation of the metal atom provided 
by oxygen is expected to be larger than that provided by hydrogen and thus pro- 
vide a better test of the quality of the RECPs. The results of this study have 
been submitted for publication, and preprints are attached [3] . 

Calculations on the platinum hydrides PtH, PtH + and PtH 2 have been car- 
ried out at the nonrelativistic (NR), perturbation theory (PT) and DHF levels. 
The DHF calculations employed the new open-shell code described above. Geome- 
tries, dipole moments, harmonic frequencies, infrared intensities and disso- 




ciation energies have been calculated. For purposes of separating spin-orbit 
effects from the non-fine-structure effects and for determining the validity 
of the first-order PT approximation, it is proposed to include in the results 
for the platinum hydrides calculations with the spin-free no-pair code of Hess 
[4] . A preliminary draft of the paper without these calculations is attached. 
These results were presented at the West Coast Theoretical Chemistry Confer- 
ence, held in May 1992 at Pacific Northwest Laboratory. 
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the effect of omission of two-electron integrals 
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The effect of omission of two-electron integrals involving basis functions for the small component of the wavefunction on the 
eigenvalue spectrum in the Dirac-Hartree-Fock problem is studied. From an analysis of the Fock matrix it is shown that omission 
of these integrals moves the negative-energy states down, not up. Their complete omission does not give rise to intruder states. 
The appearance of intruder states occurs when only some of the core integrals are omitted, due to the nature of particular contrac- 
tion schemes used for the core basis functions. Use of radi ally lo calized functions rather than atomic functions alleviates the 
intruder state problem. 


In a recent paper [ 1 ], the effect of the omission of 
two-electron integrals involving four small-compo- 
nent basis functions in Dirac-Hartree-Fock (DHF) 
calculations on some properties of PbH 4 was inves- 
tigated. The reason for the omission was a desire to 
reduce the large number of integrals involving the 
small component to a more manageable size. The 
justification for the omission was that these integrals 
only contribute to the energy at order a 4 (a is the 
fine-structure constant), and the contribution is 
mainly to the core energy. Since the effect on these 
properties was found to be small, a subsequent DHF 
study of the molecular properties of the group IV hy- 
drides [2] omitted them entirely. The negative en- 
ergy spectrum was ignored in these studies. In cur- 
rent work on PbO, the effect of successively adding 
the (SS|SS) core integrals from the Pb n— 1, 2 and 
3 shells was investigated. These integrals are few in 
number and their inclusion should recover most of 
the energetic effect of the (SS|SS) integrals. How- 
ever, in the course of this study, it was found that the 

Correspondence to: K.G. Dyall, NASA Ames Research Cen- 
ter, MS RTC 230-3, Moffett Field, CA 94035-1000, USA. 

The notation indicates that each basis function belongs to a 
given component type, L or S for large or small. 


highest negative energy state was increasing in en- 
ergy as more (SS|SS) integrals were added, and en- 
tering the region of the positive energy states. This 
observation prompted the present investigation, to 
understand why, contrary to intuition, the negative 
energy states moved up as more ( SS [ SS ) integrals 
were added, and why they crossed into the region of 
the positive energy states, above -2mc 2 * 2 . 

It is helpful in this context to consider the “neg- 
ative energy” states - those with eigenvalues nor- 
mally below - 2 me 2 - as positive kinetic energy states 
for positrons: their energies are then — (e + 2 me 2 ), 
where e is the eigenvalue. This is in line with the 
interpretation of the Dirac negative energy states 
from quantum electrodynamics (QED) [3]. The as- 
sumption that the addition of the (SS | SS) integrals 
should make the negative energy states go down is 
based on the effect of the integrals on the electron 
(“positive energy”) states, which is to increase the 
repulsion, making the electron distribution more dif- 
fuse and thereby decreasing the attraction for posi- 
trons. While this argument is correct, the reduced at- 

The electron rest energy has bee subtracted, so that the zero of 
energy corresponds to the nonrelativistic definition of a free elec- 
tron with zero kinetic energy. 
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Table 1 

Effect of addition of ( LL | SS ) and ( SS | SS ) integrals and of con- 
traction on the highest positron eigenvalue for Hg 78+ . Value is 
given in E h relative to —2 mc 2 \ basis set is tightest 14 s functions 
from 22s 1 6p 1 3d8f set of ref. [ 4 ] 


Integrals 

added 

Fully 

contracted 

Half- 

contracted 

Uncontracted 

( LL S LL ) 

-12115 

-583.16 

-580.71 

(LLjSS) 

-12009 

-572.19 

-569.75 

(SS|SS) 

-12004 

-571.06 

-568.63 


traction is not the major effect. This can be seen by 
examining the contribution the (SS | SS ) integrals 
make to the Fock matrix. It is sufficient here to con- 
sider the case where there is only one occupied spi- 
nor and one basis function for each of the large and 
small components, x L and / s . The two resultant spi- 
nors can be written 



where the superscripts indicate the electron and pos- 
itron states respectively, and a and b are the coef- 
ficients, b=(0(a ] ). The diagonal ele- 

ments of the Fock matrix can be written 

= + - ( 2 ) 

If we assume that a and b have been determined with 
the (SS | SS) integral(s) omitted, and consider the 
contribution the addition of the (SS 1 SS) integrals 
makes to the diagonal elements of F, we arrive at 

&Fee=b 4 (X S X S \X S X$) , 

&F pp =a 2 bHx s X s \X s X s ) • W 


Both contributions are positive, but the diagonal ele- 
ment for the positron state moves upwards by ( V{c 2 ) 
more than the diagonal element for the electron state. 
What has been added is an important attractive term 
for the positron states. When the ( SS | SS ) integrals 
are added, the potential is changed, increasing the 
repulsion for the electron states and increasing the 
attraction for the positron states. The same is true 
for the (LL|SS) integrals. 

As a numerical illustration of the effect discussed 
above, the Hg 78 + ion is considered. The tightest 14 
s functions were taken from Faegri’s [4] 22s 1 6p 1 3d8f 
set, and fully contracted to give the results in the first 
column of table 1, where the increase in the positron 
eigenvalue as the (SS | SS ) integrals are added is 
clearly seen. The increase for the electron eigenvalue 
is 0.5 E h , compared to 5 E h for the positron eigen- 
value. That the effect is due to increased positron- 
electron attraction is demonstrated in table 2, where 
a series of calculations with an uncontracted basis set 
for Hg [4] with increasing numbers of electrons is 
presented. The highest positron eigenvalue increases 
by 0.6 E h per electron pair as the interaction is turned 
on. Most of it comes from the (LL|SS) integrals. 
However, at this rate, the increase for neutral Hg 
would be insufficient to make the eigenvalue go above 
-2 me 2 . Note also that the highest positron eigen- 
value is independent of the number of electrons when 
only the (LL|LL) integrals are included: it is essen- 
tially experiencing only the bare nuclear potential, 
since the contributions of the (LL|LL) integrals to 
it are C r (a 4 ). Neither of these illustrations showed 
any sign of the intruder states noticed for PbO, and 
on the basis of the above arguments, the omission of 
(SS | SS) integrals should cause a downward move- 
ment of the positron eigenvalues, not upward. 


Table 2 

Effect of addition of (LL| SS) and ( SS | SS) integrals on the highest positron eigenvalue for Hg as a function of the number of electrons. 
Value is given in E h relat ive to — 2mc 2 ; basis set is taken from the 22s 1 6p 1 3d8f set of ref. [ 4 ] 


Number of 
electrons 

Eigenvalue 



Change in eigenvalue 


(LLjLL) 

■f (LLiSS) 

+ (SS|SS) 

+ (LLJSS) 

+ (SS|SS) 

total 

2 

-24.181 

-23.637 

-23.581 

0.544 

0.056 

0.600 

4 

-24.181 

-23.051 

-22.982 

1.130 

0.069 

1.199 

10 

-24.181 

-21.289 

-21.184 

2.892 

0.105 

2.997 

18 

-24.181 

-18.909 

-18.786 

5.272 

0.123 

5.395 
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Basis set contraction is one possible source of 
problems in the positron spectrum. If the core is 
heavily contracted, the small-component functions 
will not be very suitable for describing scattering 
states. The positron states closest to — 2 me 2 are com- 
posed of the most diffuse basis functions, however, 
so that contraction of the core should not greatly in- 
fluence these states. In the lower region of the pos- 
itron spectrum, contraction will decrease the num- 
ber of basis functions describing the region, giving a 
sparser discrete representation of the scattering states, 
and thus shifting the highest positron eigenvalue 
downwards. This is seen in table 1 by comparing the 
second and third columns. Calculations were per- 
formed in the half-contracted case in which the only 
(SS | SS) integrals included were for the contracted 
function. The effect on the highest positron eigen- 
value was negligible, supporting the hypothesis that 
the core contracted functions do not affect the high- 
est positron states. The issue of contraction of the 
core does not therefore seem to address the problem 
that for PbO there were eigenvalues which were above 

— 2 me 2 by hundreds of hartrees when the core 
( SS | SS ) integrals were added. 

In itself, the occurrence of positron states above 

— 2 me 2 is not unphysical. It indicates bound posi- 
trons, just as electron states below zero energy in- 
dicate bound electrons. That atoms or molecules 
should be able to bind positrons should not be sur- 
prising, considering their ability to bind protons. In 
calculations on the F“ ion, states at about 0.1 E h 
above —2 me 2 were observed - a much more realistic 
value for a positron affinity. A similar observation 
was made by Visscher et al. [5] for EuO^ - , where, 
as for F - , the long-range potential is attractive for 
positrons. The size of these “positron affinities” for 
PbO, 400 £ h , is much too great to be a physical effect. 

States significantly above —2 me 2 have also been 
noted by Visscher et al. [5] when the Aitken ex- 
trapolation procedure was used in the SCF process 
to accelerate convergence, but at convergence the 
states were no longer above — Imc 1 . Several differ- 
ent methods for converging the PbO calculations, in- 
cluding the (SS | SS ) integrals for the Pb Is, 2s and 
2p orbitals only, were tried: level shifting, DIIS [6] 
extrapolation on the off-diagonal blocks of the Fock 
matrix or on the differences in the coefficients be- 
tween iterations, damping on the density matrix or 


on the coefficients, and various combinations of 
these. None had any effect on the intruder state. 

Because of the resources needed for complete cal- 
culations on PbO, further tests were carried out on 
SnO. Three different calculations were performed. 
In the first, the (LL | SS) and (SS | SS) integrals in- 
volving the small-component valence basis functions 
were omitted, in the second only the (SS | SS ) inte- 
grals over the valence basis functions were omitted, 
and in the third all integrals were included. In the 
first calculation, the highest positron eigenvalue was 
at 1 53.56 E h above — 2mc 2 , and was not the only one 
above —2 me 2 . The two highest positron states had 
large coefficients from the Sn core s functions, and 
were the highest positron states composed of these 
functions. Similarly, the PbO intruder had large 
coefficients of the Pb core s functions, and was the 
highest positron state composed of these functions. 
In the second calculation on SnO, which used the 
vectors from the first calculation as a starting guess, 
the highest positron state was 0. 1 06 E h below - Imc 2 , 
and it had only small contributions from the Sn core 
s functions. In the third it rose slightly, to 0.081 E h 
below — 2mc 2 , in accordance with the behaviour ex- 
pected from the arguments presented above. 

The evidence to this point suggests that the com- 
plete omission of a class of integrals - such as 
(SS | SS) — does not produce intruder states, but that 
addition of only some of the integrals of a given class 
may do so. In the calculations on PbO, the first stage 
included the (LL|LL) and (LL|SS) integrals. The 
addition of the (SS | SS ) integrals for the Pb n= 1 and 
2 shells produced upward shifts in elements of the 
Fock matrix in the atomic spinor basis of the order 
of 10 £ h , a not unreasonable shift. This shift was 
magnified to 90 E h for the highest core s-like posi- 
tron state in the molecular spinor basis due to the 
large coefficients in the trial vectors, which were taken 
from the first stage of the calculation. Clearly, the 
( SS | SS ) integrals which were added make a signif- 
icant contribution to the Fock matrix and cannot be 
omitted at any stage of the calculation. The absence 
of these integrals in the first stage of the calculation 
produces a distortion in the positron eigenvectors 
which, when the integrals are put back in, gives rise 
to the intruder states. The same appears to be true 
of the valence (LL|SS) integrals in the calculations 
on SnO. While the omission of the valence (SS| SS) 
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integrals for SnO appears to have had very little ef- 
fect on the positron spectrum, it remains to be seen 
whether this will be true for heavier systems such as 
PbO. It should be pointed out that for PbO, only 
some of the core (SS | SS) integrals were put back into 
the calculation. If the remaining core (SS j SS ) in- 
tegrals were included, the intruder could disappear 
as it did for SnO. It is clear also that the procedure 
in which the classes of integrals, (LL|LL), (LL|SS) 
and (SS|SS), are added in stages has some role to 
play in the appearance of intruder states. A better, 
but more time-consuming, procedure would be to use 
all the integrals and build up the electron occupation 
from the core. 

The relation of the appearance of intruder states 
to the large size of some of the coefficients in the mo- 
lecular spinor expansion re-opens the question of core 
contraction. It may be that the inflexibility of the 
general contraction in the core does not permit a suf- 
ficiently good description of the corresponding pos- 
itron states to avoid the appearance of intruders; it 
may be that all of the core integrals are necessary to 
ensure that the contracted functions do give a good 
description of the positron states, or it may be sim- 
ply that certain classes of integrals cannot be omitted 
regardless of the nature of the contraction. Addition 
of extra functions for the small component in order 
to give it a better description is to be avoided for sev- 
eral reasons. The first is that the additional functions 
will add significantly to the cost of any calculations 
- not a small consideration, given the large number 
of integrals required. The second is that generally 
contracted core functions taken from atomic calcu- 
lation should provide a good description of the mo- 
lecular core, and thus any additional functions will 
only go into the description of the positron states, 
which are not of interest for most molecular appli- 
cations. The third is the matching of functions for 
small and large components and their one-to-one 
correspondence with electron and positron states, re- 
quired by physical considerations [7]. 

An alternative to the current general contraction 
scheme, which uses the atomic SCF functions, is to 
transform the core basis functions to a form in which 
they are readily localized - “nodeless” functions. 
While the same space would be covered by these 
functions, the effect of inclusion of a small subset of 
the (SS | SS ) integrals from these functions in cal- 


culations such as were done for PbO should be mark- 
edly different. The inner tail of the core functions is 
the part from which the largest (SS|SS) integral 
contribution comes: in the atomic contraction 
scheme, each function has its own inner tail, whereas 
in the radially localized contraction scheme, there is 
only one function for the inner tail. Thus, adding 
(SS | SS) integrals for only a few core functions in 
the atomic scheme creates a very uneven distribu- 
tion of the inner tail contribution, but in the radially 
localized scheme, most of the inner tail contribution 
is obtained, and it is more evenly distributed. 

To test this hypothesis, radially localized functions 
(RLFs) were generated from the atomic functions in 
a Gaussian basis as follows. Each RLF for a given n 
quantum number was expressed as a linear combi- 
nation of the atomic function and the RLFs with 
lower n value: 

^ oc =0„- "icrfdr. ( 4 ) 

k = 1 

Each function is expressed as a linear combination 
of Gaussian functions g M , 

/V' .v 

<Pk c = X b kM g M , <p k = I a k „g M , (5) 

M= 1 

so that the RLF can be written 

<t>Y~ X [a,,,,- X Cn k b k Ag» ■ (6) 

A convenient way of performing the radial locali- 
zation is to minimize the coefficients of the Gaus- 
sian expansion of the RLF over some range of ex- 
ponents. We choose the function to be minimized as 
the weighted sum of the squared coefficients, 



with weights w,,=JT», where ^ is the exponent of 
Gaussian ji , and the cutoff value for the sum, n c , is 
chosen to correspond to the ( /? — 1 ) th radial node, 
both the atomic functions and the radially localized 
functions have been used in a series calculations on 
PbO at r=3.6 a 0 * The first stage of the calculations 
included only the (LL|LL) and (LL|SS) integrals. 
At each succeeding stage, ( SS | SS ) integrals were 
added for the next highest n quantum number on Pb, 
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Table 3 

Effect of contraction scheme on the highest positron eigenvalue and the total energy and energy change as successive shells of (SS)SS) 
integrals are added to calculations on PbO at 3.6 a Q . Eigenvalue is given relative to —2 me 2 . First row has results without (SS|SS) 
integrals. All energies are given in E h 


(SS|SS) 

integrals 

added 

Atomic SCF contraction 


Radially localized contraction 


eigenvalue 

E 

Mot 

ATtot 

eigenvalue 

ftot 

AE'ioi 


-0.1583 

-20989.59304 


-0.1583 

-20989.59304 


n= 1 

-0.1579 

-20989.02265 

+ 0.57039 

-0.1583 

-20988.68970 

+ 0.90334 

n = 2 

4- 128.06 

-20986.72719 

+ 2.29546 

-0.1578 

-20985.36697 

+ 3.32273 

n = 3 

+ 468.55 

-20984.96186 

+ 1.76533 

-0.1576 

-20984.39280 

+ 0.97417 

n~4 

+ 619.94 

-20984.85034 

+ 0.11152 

-0.1569 

-20985.22919 

-0.83639 

n = 5 

-0.1580 

-20986.74701 

-1.89037 

-0.1580 

-20986.74071 

-1.51152 


starting with Pb n— 1, and adding the O n= 1 and 2 
with the Pb n = 4 and 5 integrals. The results are pre- 
sented in table 3, labelled with the Pb n quantum 
number. Three features emerge from this study, con- 
firming the hypothesis made above. The first is that 
there are no intruder states in the series of calcula- 
tions with the RLFs. The highest positron eigenvalue 
displays a small upward shift as the (SS|SS) inte- 
grals are successively added. The second is that a 
larger energy contribution is obtained from the lower 
n (SS | SS) integrals using RLFs than using the atomic 
functions. Both of these indicate that the contribu- 
tions for all the core functions are being added in or- 
der of importance in the RLF series of calculations. 
The third is that at n- 5, the two calculations give 
the same results, as expected, and no intruder states 
are in evidence. 

The conclusion of this investigation is that the ap- 
pearance of the intruder states was related not to 
general contraction per se, or to the SCF conver- 
gence method, but to the particular general contrac- 
tion scheme used. To ensure that they do not appear, 
either all of the core (SS | SS) integrals have to be 
included in the calculation, or if this is not feasible, 
a radially localized contraction scheme should be 
employed to ensure a more even distribution of the 
( SS | SS ) contributions to the Fock matrix. The cal- 


culations on SnO indicate that it is possible to per- 
form meaningful calculations which omit the va- 
lence ( SS | SS ) integrals. However, care must always 
be exercised in the omission of classes of integrals, 
and there may be cases in which it is not possible to 
omit subsets of integrals involving the small com- 
ponent. This would restrict the range of problems 
which could be tackled with the DHF method using 
conventional SCF procedures and current storage 
availability. 

The author was supported by NASA grant No. 
NCC2-552, and would like to thank Dr. Winifred 
Huo and Dr. Peter Taylor for helpful discussions. 
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Abstract 

The 2 S + _ se p ar ation in KO is investigated using large basis sets and high 
levels of correlation treatment. Relativistic effects are included at the Dirac-Fock 
level and reduce the separation only slightly. The basis set superposition error is 
considered in detail. On the basis of these calculations, our best estimate places 
the 2 H 3/2 state about 200 cm -1 above the ground 2 £ + state in agreement with our 
previous estimate. 

I. Introduction 

In 1986 we computed the 2 £~^" — “II separation in KO to be 238 cm with 
the 2 E+ state being the lower of the two [1]. This separation was computed at 
the singles and doubles configuration interaction (SDCI) level using a large Slater 
basis set. While not definitive, the separation was thought to be accurate because 
of the small difference between the SDCI and self-consistent-field (SCF) values. 
However, when 15 electrons are correlated, an SDCI treatment can underestimate 
the differential correlation effects. Approaches such as the modified coupled-pair 
functional (MCPF) method [2], or the coupled cluster singles and doubles method 
with a perturbational estimate of the triple excitations [3] [CCSD(T)], account for 
the effect of higher than double excitations and can therefore yield superior results 
to the SDCI approach. 

Experiment all y, the ground state has not been determined and there is conflict- 
ing evidence. The failure to observe an electron spin resonance (ESR) spectrum [4] 
is consistent with a 2 H ground state, while the magnetic deflection experiments [5] 
suggest a 2 £ + ground state. Photoelectron detachment experiments [6] on KO 

f Mailing address: NASA Ames Research Center, Moffett Field, CA 94035-1000 
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might offer a mechanism to directly observe both states, and hence determine the 
ground state and the separation; such experiments [6] have recently been performed 
for NaO - and preliminary results suggest that they are possible for all of the alkali 
oxides systems. In addition, microwave experiments [7] have recently been per- 
formed on NaO and, with a pure-precession hypothesis, an estimate of the 2 E~ — 2 II 
separation was obtained. An analogous experiment should be possible for KO. 

Given the recent experimental activity on the alkali oxides, it seems timely to 
reinvestigate the 2 £ + - 2 II separation in KO at higher levels of theory. We consider 
the effect of one-particle basis set and basis set superposition error (BSSE) in more 
detail than in our previous work. In addition, the separation between the 2 E + and 
2 II states is small enough that relativistic effects could affect the ordering of the 
states. These were not included in the previous work and are addressed here. 

II. Methods 

Two basis sets are used in this work. The small oxygen basis set is the 
(135 9 p 6d 4/)/[5s 5 p (2 + 1 )d If] atomic natural orbital [8] (ANO) set [9] de- 
termined from the average natural orbitals of 0 and 0 . (The notation n + m in- 
dicates n ANO orbitals are employed and m primitive functions are uncontracted). 
In the large oxygen basis set only two s and one p ANOs are used. The outer 
five 3 and five p primitives and the d and / polarization functions are uncon- 
tracted. This is augmented with two uncontracted g functions and a diffuse s 
function with an orbital exponent of 0.076666. Thus the final basis set is of the 
form (145 9p 6d 4/ 2g)j{{2 + 6)s (1 + 5 )p 6d 4/ 2 g\. 

The small K basis is derived from the (245 15 p) basis [10] augmented with an 
5 function optimized for K~ and three p functions optimized for the 2 P state of K. 
The 5 and p basis sets are contracted to (3+6) and (2+7), respectively, where the 
three s and two p contraction coefficients are taken from the SCF wave function. 
The basis is further augmented with six even-tempered (a n =2.5 7l ao for n= 0,fe) d 
functions with ao =0.0455 and four even-tempered / functions with olq =0.12949. 
The small K basis set is of the form (255 18p 6 d 4/)/[(3 + 6)5 (2 + 7)p 6 d 4/]. In 
the large K basis set, the 5 and p functions are contracted (2+12) and (2+9) and 
two g functions (a=0.768 and 0.307) are added, thus yielding a basis set of the form 
(25 5 18p 6d 4/ 2g)/[(2 + 12)5 (2 + 9)p 6 d 4/ 2 g]. Only the pure spherical harmonic 
components of the basis functions are used. 
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The orbitals are optimized using the SCF approach. We have treated the 2 II 
state using symmetry and equivalence restrictions (i.e., a 7 t£' 5 7t 2 ' 5 occupation ) or 
by using a symmetry broken solution (i.e., a 7 r 2 7 r 2 occupation). Electron correlation 
is added using the MCPF or CCSD(T) approaches. Fifteen electrons are correlated, 
corresponding to the 0 2s and 2 p and the K 3s, 3p, and 4s electrons. 

Relativistic effects are estimated as the difference between the SCF and Dirac- 
Fock (DF) results. Because correlation was not included, it was possible to reduce 
the basis set size without affecting the accuracy of the results. The full s and p 
primitive basis sets for K and 0 were contracted to [(3 + 6)s (2 + 7 )p] and [(2 + 
3)s (1 + 4)p] using the orbitals from the SCF wave functions for K + and 0~. For 
polarization the 6 d set and the outer three of the 4/ functions on K were used, and 
on 0 the outer 3 primitive d functions and the / function with the largest coefficient 
from the ANO polarization set were used. The SCF 2 X) + — 2 II separation in this 
basis set was identical to that found using the small ANO set. 

The SCF/MCPF calculations were performed using the SEWARD [11]- 
SWEDEN [12] program system. The CCSD(T) calculations were performed using 
the program developed by Scuseria [13]. The atomic DF calculations were performed 
using GRASP [14] while the molecular calculations were performed using the code 
developed by Dyall [15]. All calculation were performed on the NASA Ames Central 
Computer Facility CRAY Y-MP or Computational Chemistry Branch IBM RISC 
System/6000 computers. 

III. Results and Discussion 

The separation at the SCF level is only slightly larger than in our previous 
work [1], with the big basis set result essentially reaching our previous estimate of 
250 cm -1 for the numerical Hartree-Fock separation. The SDCI separation in the 
small basis set is only 6 cm -1 smaller than that reported in our previous study. 
With a symmetry broken SCF solution, the 2 13 state is predicted to be below the 
2 E + state. When correlation is added at the SDCI level, the 2 E" r state is predicted 
to be lower. However, the SDCI based on the symmetry broken SCF orbitals yields 
a significantly s mall er separation. This suggests that the SDCI approach is unable 
to overcome the bias in the SCF wave function. At the MCPF level, the effect of 
symmetry breaking is smaller, as expected, since the MCPF approach includes the 
effect of higher excitations. At this level, the symmetry broken solution actually 




leads to a larger separation than using the SCF with symmetry and equivalence 
X restrictions. Both MCPF separations are about 100 cm 1 smaller than the SDCI 

result with symmetry and equivalence restrictions. The CCSD separation is the 
same as the MCPF. The effect of triples is to further increase the separation by 
47 cm -1 , which is almost an order of magnitude smaller than the difference be- 
tween the SCF and CCSD (368 cm -1 ). This strongly suggests that even higher 
levels of theory will not significantly increase the correlation contribution to the 
separation. Expanding the basis set results in a small increase in both SCF and 
MCPF separations. 

In Table II we report the BSSE at the MCPF level. We compute the BSSE for 
K + with the 0 basis set, K + (0), and O - with the K basis set, because both states 
of the molecule at r e are best described as K + 0 - . K + is a closed shell; therefore the 
calculation of the BSSE is straightforward, and we consider K' r (0) at the two rele- 
vant internuclear separations. The 2p 5 occupation of 0 is oriented 2p<7"2p7T 3 and 
2p(7 1 2p7r 4 for the 2 II and 2 E + states, respectively. As shown in previous work [16], 
the BSSE can differ significantly for the two occupations. Therefore we use the 
BSSE for the 2p<7 2 2p7r 3 occupation to correct the 2 II state and the BSSE for the 
2p<x 3 2p7r 4 occupation to correct the 2 S"*’ state. 

For K + the BSSE is small for both basis sets. It is slightly larger for the 2 E+ 
state than the 2 II state because of the shorter bond length. In contrast, the 0“ 
BSSE is much larger for the smaller basis set. In spite of the longer bond length, 
the 0 - (K) BSSE is larger for the 2 II state, because a pair of oxygen 2p<r electrons 
points at the ghost basis set compared with only one electron in the 2 E + state. 
For both basis sets the total BSSE is very similar for the 2 E + and 2 II states and 
therefore does not have an important differential effect on the separation. Given 
the small total BSSE for the big basis set, this conclusion seems definitive in spite 
of any uncertainties in how to best compute the BSSE. 

One effect not considered to this point is relativity. The spin-orbit splitting in 
the 0" ion was calculated to be 193 cm -1 at the Dirac-Coulomb level. Inclusion 
of the Breit interaction decreased the splitting to 170 cm 3 . This value is expected 
to be accurate as an analogous treatment of O agrees with experiment to within 
2 cm -1 . Because the bonding in KO is best described as K"*"0 , this splitting in 0 
can be used in a simple first-order model [17] to estimate the splitting in KO. On 
this basis, the 2 E^ 2 state is unshifted from its nonrelativistic position, the 2 TL 1 / 2 
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state is shifted up by 57 cm -1 and the 2 n 3 / 2 state is shifted down by 57 cm" 1 . 

Mixing of K 45 character from covalent contributions to the wave functions 
and spin-orbit mixing of the 2 H \/2 and 2 E^ 2 states could invalidate the first-order 
model. Therefore we study the 2 n 3 /2 and 2 2^ states at the Dirac-Fock level. 
On the basis of the nonrelativistic calculations we use a bond length of 4.4 Go for 
the 2 n 3/2 state and 4.13 Go for the 2 S^ 2 state. At this level, the 2 n 3 state 
is lowered by 36 cm -1 and the 2 S^ 2 state is lowered by 24 cm -1 , relative to the 
nonrelativistic states. Thus the differential effect of relativity on the state separation 
is only 12 cm" 1 , or about 1/5 of the estimate derived from 0”. Clearly, however, 
the differential effect of relativity is sufficiently small that it will not reverse the 
order of the states computed at the nonrelativistic level. 

Our best separation is obtained by correcting the CCSD(T) separation for: 1) 
MCPF basis set improvement (+36 cm" 1 ), 2) BSSE (+8 cm" 1 ) and 3) relativistic 
effects (-12 cm" 1 ). This yields a separation of 214 cm" 1 , with the 2 E + state 
being lower. This is very similar to our previous result [1], While it is difficult 
to be definitive for a separation this small, the calculations strongly suggest that 
the ground state is 2 E + . Because improving the basis set increases the separation, 
as does including higher levels of correlation treatment, and relativistic effects are 
very small, it is difficult to imagine what effect would reverse the order of these two 
states. 

IV. Conclusions 

The separation at our highest level of theory, the CCSD(T) approach, is 
182 cm" 1 . Including the correction for basis set incompleteness from the MCPF 
level increases the separation to 218 c m" 1 , which is very similar to our previous 
result [1]. The differential effect of BSSE is 8 cm" 1 , increasing the separation. Rel- 
ativistic effects reduce the separation by only 12 cm" 1 . Thus we conclude that the 
ground state is 2 E + . 
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Table II. Summary of the BSSE at the MCPF level. 


BSSE(cm *) 

Small basis set 

2 S + 

K+(0) at r=4.13 32.1 

0~(K) at r=4.13 348.4 

Total 380.5 

2 n 

K+(0) at r=4.40 20.9 

0~(K) at r=4.40 364.5 

Total 385.4 

Big basis set 

2 S+ 

K+(0) at r=4.13 32.5 

0~(K) at r=4.13 88.5 

Total 121.0 

2 n 

K + (0) at r=4.40 25.8 

0“(K) at r=4.40 103.4 

Total 129.2 
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All-electron molecular Dirac-Hartree-Fock calculations: 
properties of the Group IV monoxides GeO, SnO and PbO. 
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California 94303, TJ.S.A. 


Abstract 

Dirac-Hartree-Fock calculations have been carried out on the ground states 
of the group IV monoxides GeO, SnO and PbO. Geometries, dipole moments and 
infrared data are presented. For comparison, nonrelativistic, first-order pertur- 
bation and relativistic effective core potential calculations have also been carried 
out. Where appropriate the results are compared with the experimental data and 
previous calculations. Spin-orbit effects are of great importance for PbO, where 
first-order perturbation theory including only the mass-velocity and Darwin terms 
is inadequate to predict the relativistic corrections to the properties. The rela- 
tivistic effective core potential results show a larger deviation from the all-electron 
values than for the hydrides, and confirm the conclusions drawn on the basis of the 
hydride calculations. 


* Mailing address: NASA Ames Research Center, MS RTC 230-3, Moffett Field, 
CA 94035-1000, U.S.A. 



I. INTRODUCTION 


Relativistic effects are known to be important in the chemistry of the heavy 
elements, but how important is an issue which can only be decided on the ba- 
sis of accurate calculations. The demands of calculations on molecules containing 
heavy elements including relativistic effects are such that until recently, various ap- 
proximate methods had to be employed to reduce the size of the calculations. All 
such methods are based in some way on the Dirac equation, which is in turn an 
approximation to equations derived from quantum electrodynamics [1], 

The most common approximation is the use of a relativistic effective core 
potential (RECP) [2,3] in which the Dirac equation itself or an approximation to it 
such as the Cowan-GrifRn equation [4] is used to generate valence pseudo-orbitals 
from which the effective potential is constructed. Usually, the spin-dependent terms 
are averaged out to give an RECP which can be used with standard nonrelativistic 
codes. Several sets of RECPs have been published [5—8,9]. 

Some other approaches, which do not depend on the frozen-core approxima- 
tion as do the RECPs, are first-order perturbation theory (PT) with the spin-free 
terms — the mass-velocity and Darwin (MVD) terms — in the perturbation Hamil- 
tonian [4,10], and the spin-free no-pair method of Hess et al. [11,12]. These are 
based on a transformation of the Dirac equation to eliminate the small component 
to a given order and truncation of the r esulta nt expression to obtain a spin-free 
Hamiltonian. The method of Hess et al. provides a Hamiltonian which can be used 
in variational calculations, whereas the MVD operator is strictly a perturbation 
operator. 

The past few years has seen the development of some all-electron Dirac- 
Hartree-Fock (DHF) codes [13-15]. Methods for inclusion of electron correlation 
based on these codes are only now being developed [16]. Although correlation effects 
are important in obtaining quantitative predictions of molecular properties, it is 
important to calibrate the various approximations to the Dirac equation at the self- 
consistent field (SCF) level, because if the approximate methods are inaccurate at 


1 



the SCF level, they can only provide correlated results in agreement with experiment 
by accident. 

This paper is the third of a series exami ning relativistic effects on proper- 
ties of small molecules containing Group IV elements, and providing calibration of 
RECP and PT methods. The first and second papers [17,18], hereafter referred to 
as I and II, were concerned with the hydrides. Despite the shortness of the bond 
lengths, hydrogen offers only a small perturbation to the central atom. Introducing 
a heavier atom should give a better test of the quality of the approximate methods. 
In this paper the properties of the monoxides of Ge, Sn and Pb are examined us- 
ing nonrelativistic (NR) SCF, PT including only the MVD terms, RECPs, and the 
DHF method. 

The monoxides are well-known experimentally [19]. Previous calculations 
have been performed on PbO by Basch et al. [20] at the SCF and multi-configuration 
self-consistent field (MCSCF) levels of theory using RECPs, by Schwenzer et al. [21] 
at the NR SCF level, by Datta et al. [22] at the SCF level using both relativistic 
and nonrelativistic ECPs, and by Balasubramanian and Pitzer [23] at the singles 
and doubles configuration interaction (SDCI) level with an RECP for lead and 
including spin-orbit interaction at the Cl stage. Balasubramanian and Pitzer have 
performed similar calculations SnO [24]. Igel-Mann et al. [25] have studied SnO and 
other molecules extensively, at the SCF, complete active space self-consistent field 
(CASSCF) and SDCI levels using RECPs with a core polarization potential and 
several basis sets. Bouteiller et al. [8] have performed calculations on GeO (among 
other molecules) at the SCF level, both all-electron and with an ECP. Comparisons 
are made with these calculations where appropriate. 


II. COMPUTATIONAL DETAILS 

The primitive basis sets used for the group IV elements were the same as 
used in I and II. The exponents were energy-optimized in nonrelativistic SCF cal- 
culations and are of approximately valence double-zeta quality. The d basis set was 
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supplemented with 2 extra functions to describe valence polarization. The primitive 
basis set for oxygen was the van Duijneveldt [26] lls6p basis supplemented with a 
diffuse s and p function to help describe the negative-ion character in the molecules, 
and 2 d polarization functions taken from Dunning’s [27] pVTZ basis. Contraction 
coefficients for all basis sets were determined from atomic SCF calculations using a 
general contraction scheme. The relativistic contraction coefficients were obtained 
from an adaption of GRASP [28], as described in I. For the Group IV elements, all 
functions up to the (n — l)d shell were kept in the core, and the outermost three 
5 , p and d primitive functions were uncontracted to form the valence basis. The 0 
basis was contracted to 5s4p2d, with three uncontracted s and p functions in the 
valence basis, and the inner tails of the 2s and 2 p functions in the core. 

Four sets of RECPs were used for Ge, Sn and Pb: those of Hay and Wadt 
[5] (hereafter referred to as HW), those of Stevens, Krauss, Basch and Jasien [6] 
(hereafter referred to as SKBJ), and both the full-core and the semi-core potentials 
of Ref. 7 (hereafter referred to as CER and CER+d respectively). Only the spin- 
averaged potentials (averaged relativistic effective potentials, AREPs) from ref. 7 
are used in this work. Both the valence sp basis sets from the all-electron calcu- 
lations and the sp basis sets supplied with the RECPs, contracted to 3$3p, were 
used in the RECP calculations, supplemented by the valence d functions from the 
all-electron basis set. The d orbital supplied with the CER+d basis was left fully 
contracted. 

The equilibrium geometry of the molecules was determined from a quartic 
fit to 5 or 6 points around r c . The energy at the predicted r e was added to the fit 
to determine the force constants. A quartic fit to the dipole moments was also used 
to determine the dipole derivatives at r e . The program INTDER [29] was used to 
obtain the harmonic frequencies and infrared intensities. The isotopes used for the 
frequencies were 16 0, 74 Ge, 120 Sn and 208 Pb. All properties are reported at the 
predicted r e value for each method. The value of the speed of light in atomic units 
was taken to be 137.03604. 

The MOLECULE/SWEDEN [30] package was used to obtain the NR SCF, 
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PT and RECP results. DHF results were obtained with the program described 
previously [15,17]. To limit the size of the DHF calculations, it has been found useful 
to discard all of the integrals involving the small component valence basis functions 
[17,31]. This has a negligible effect on the properties even for Pb compounds, 
because the terms omitted contribute the energy only at 0(a 4 ). In the present 
calculations these (SS|SS) integrals were discarded for PbO. All calculations were 
performed on the Computational Chemistry Branch Convex C-210 and the Central 
Computing Facility CRAY Y-MP/864 computers at NASA Ames Research Center. 


III. RESULTS AND DISCUSSION 

The calculated bond lengths r e are presented in Table I, along with the 
experimental values [19], and the relativistic corrections to the bond lengths A rel r e 
predicted by the PT and DHF calculations. The dipole moment data are presented 
in Table II, and the harmonic frequencies and infrared (IR) intensities in Table III. 
The first part of the discussion deals with the all-electron data and the properties 
of the molecules in general, in the second part the RECP data are discussed and 
compared with the all-electron data, and in the third, comparisons are made with 
other calculations. 

A. All-electron results. 

An understanding of the differences in properties of the group IV monoxides 
is aided by consideration of the relative energies of the atomic orbitals. The orbital 
energies for Ge, Sn, Pb and O obtained from configuration average SCF calculations 
are given in Table IV. The O 2s orbital is close to the d orbital on the meted, and 
the metal s orbital is close to the 0 2 p orbital. It is therefore expected that the 
O 2s orbital will not participate much in the bonding. The table also shows that 
Sn has the least bound valence orbitals when relativistic effects are included (the 
6p 3 / 2 is unoccupied in the ground configuration of Pb): experimentally, Sn has the 
smallest ionization potential (IP) in the group. Sn may thus be expected to give 
the most ionic oxide. This conclusion would not be borne out by nonrelativistic 
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calculations, for which the IP decreases monotonically down the group. 

The relativistic bond length corrections are smaller than those of the hy- 
drides. Most of the bond length contraction resulting from the spin-free terms in 
the relativistic Hamiltonian is cancelled by the spin-orbit interaction, which causes 
a partial promotion of a electrons into tt orbitals that are usually of antibonding 
character. Comparison of the DHF bond lengths with experiment shows a uni- 
form underestimate of 0 . 03 A, partly due to basis set effects and partly to electron 
correlation. 

The dipole moments mostly show very little change from relativistic effects. 
There is a small increase for GeO and SnO, and the spin-free terms give a small 
decrease for PbO. The dipole moment of PbO is reduced significantly at the DHF 
level, because of significant changes in the valence molecular spinors (MSs) from 
their nonrelativistic counterparts. PT is unable to describe these changes, and 
thus does not give an accurate value for The DHF atomic charges obtained 
by dividing fi e by r € show SnO more ionic than GeO by 0.07e, and PbO less ionic 
than SnO by O.Ole. Without spin-orbit interaction, PbO is more ionic than SnO 
by 0.02e. It should be noted that the dipole moments reported here include effects 
from both geometric and electronic structural changes due to relativity. If the dipole 
moments at a fixed geometry are considered, the changes are larger. For PbO at 
3.6ao, for instance, PT gives a relativistic correction of 0.27 D, compared to the 
DHF correction of —0.13 D. The spin-orbit effect at 3.6ao is then —0.40 D, rather 
than the —0.20 D obtained by comparing results at the respective r e values. 

It is for the IR harmonic frequencies and intensities that the effect of spin- 
orbit interaction is most dramatic, resulting in a 10% reduction in the frequency for 
PbO, and more than halving the intensity. The changes in the valence MSs bring 
about a large decrease in the dipole derivative. PT obtains only 7% of the frequency 
reduction, and predicts an increased intensity. The frequency reduction is smaller 
for SnO but PT still obtains only a fraction of it, whereas for GeO spin-orbit effects 
are small enough that PT gives a satisfactory result. For the intensities of GeO and 
SnO, spin-orbit interaction apparently has no effect. 
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It is clear that spin-orbit interaction has a strong influence on the properties 
of PbO and by inference on the bonding. The nature of the bonding in heavy p- 
block diatomic molecules is strongly influenced by spin-orbit interaction and must 
be discussed in terms of w-u? coupling. Several authors have discussed aspects of the 
bonding for the p - block elements [32-34]; the discussion is presented and extended 
here. Both p 1 j 2 and p 3 / 2 atomic spinors may contribute to w = 1/2 molecular 
spinors. The form of these atomic spinors (with rrij = = 1 / 2 ) is 


1 1 

f -pa \ 

1 / 

f s/2pa\ 

p,/2 - 73 ' 

^y/2jm) ’ 

73/2 " 75 1 

\ PTT / 


(i) 


An a? = 1/2 molecular spinor may be written 

/2 = C 1 Pi /2 "b C 2 Pi /2 "b c 3 P3 /2 "b c 4 P3 J2 

_ / pa A (\Z2c 3 - Cl) + pcr B (v / 2c 4 - c 2 ) 
Vp7T >1 (v/2ci + C 3 ) + p7r B (v / 2c 2 + c 4 ) 


The bonding character of the spinors will be determined by the coefficients. 


For homonuclear diatomics, ci = ±c 2 and C 3 = ±c 4 for the u and g combi- 
nations, giving 


_ (P<r g {V2c 3 - Ci)\ 


_ / p<t u (\/2c 3 - Ci)\ 

^ / 2 u "U 7 r u (^ci+c 3 ); 


( 3 ) 


The molecular symmetry thus imposes on the spinors a mixture of bonding and 
antibonding character, which is determined by the ratio of the coefficients c\ and 
C 3 . The A -s coupling limit is obtained with c\ = y/2c$ for a pure 7 r/? spinor or 
c 3 = — \/2ci for a pure aa spinor. This requires a promotion from the P\j 2 to 
the P 3/2 atomic spinor, which for light elements is negligible, but for heavy ele- 
ments can be considerable, resulting in weaker bonds. In TI 2 for instance, the 
promotion cost is sufficient to render the lowest state bound by only O.OleV at 
the 2-configuration SCF level [33], and the ground state is the 0“ state from the 
configuration (l/2 p ) 1 (l/2 ti ) 1 , bound by only 0.04eV. 
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For heteronuclear diatomics, the coefficients are not restricted by the gju 
symmetry as they are for homonuclear diatomics, and hence the u— bonding/ 7 r— 
antibonding and <r-antibonding/7r-bonding combinations need not arise. The re- 
strictions will instead be on the ratios of cj to C2 and C3 to C4, determined by 
the magnitude of the spin-orbit splitting of the p shell on each centre. For small 
spin-orbit splittings the \-s coupling limit is obtained with coefficients as for the 
homonuclear case. In the extreme of large spin-orbit splitting on centre A, it is 
possible to obtain pure bonding spinors with C3 = C2 = 0: 


^ 1/2 = 


/ -per A + \/ 2 ; p(TB 

V V2pnA + P^B 


( 4 ) 


In this spinor, the cr component is skewed towards centre B and the 7T component 
towards centre A. This kind of combination of atomic spinors would be favoured in 
SCF energy optimization. The opposite skewing is obtained with cj = C4 = 0, but is 
not favoured for large spin-orbit splitting on centre A due to the promotion cost. 


The DHF Mulliken gross populations presented in Table V and the integrated 
spin densities presented in Table VI illustrate well the trend from small to large 
spin-orbit interaction. For GeO with small spin-orbit interaction, the pi/ 2 : P3/2 
ratios are close to the 1:2 ratio required for a <r orbital with a spin in the 12e!/ 2 
and 14e!/ 2 spinors and 2:1 for a 7r orbital with (3 spin in the 13e!/ 2 spinor. These 
spinors are spin-pure to better than 1 part in 1000. For SnO, in which the spin-orbit 
interaction is of moderate strength, the ratios are close to 1:1 on the Sn atom and 
approximately 3:1 and 1:4 on the O atom for the 18ej/ 2 and 19ej^ 2 spinors. The 
17e 1/2 spinor still has the LS ratio. PbO has large spin-orbit interaction, and the 
25ei/2 and 26ej/2 spinors have either pi/ 2 011 O with p3/ 2 on Pb or vice versa. Even 
the 24ej/2 spinor has a 1:1 p\/i : P3/2 ratio on the O atom. The 2he- l /2 and 26ej/2 
spinors have almost equal a and /3 spin densities. The smaller 6p 3 / 2 population in 
the 25ej/ 2 spinor compared to the 6pi/ 2 population in the 26e!/ 2 spinor reflects the 
promotion cost from the 6pi / 2 to the 6p3 / 2 atomic spinor. 

A further possible source of spin-orbit effects is the underlying d shell, which 
has a larger spin-orbit splitting than the valence p shell, and which is close in 
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energy to the 0 2a shell. Despite this near-degeneracy and the resultant mixing of 
the atomic spinors that induces non-statistical 5d 3 / 2 :5d5/ 2 ratios in the molecular 
spinors, the sum over the three relevant cv = 1/2 spinors gives the statistical ratio, 
and there is no net spin-orbit effect on the molecule from the 5 d shell. 

The nonrelativistic Mulliken gross populations are given in Table VII. Com- 
parison with the DHF populations shows a slight relativistic decrease in the atomic 
charges, a small decrease in the metal s population and increase in the 0 s popula- 
tion, with corresponding changes in the p populations. At first sight, the decrease 
in the metal s population is counter to the expected effect of relativity. The contri- 
butions from the individual orbitals (or spinors) gives some insight into this effect. 
The lower of the valence cr orbitals is principally a bonding combination between 
the metal s orbital and the 0 2 pa orbital, and the higher orbital an antibonding 
combination with an sp hybrid on the metal polarized away from the 0 atom. Due 
to the relativistic stabilization of the metal s orbital, its contribution to the lower <r 
orbital increases. Consequently, its contribution to the higher a orbital decreases, 
with the overall effect of decreasing the total metal s population. The effect is more 
pronounced for PbO than the lighter oxides. The spin-orbit stabilization of the 
6p 1/2 spinor in Pb assists in the transfer of charge from the 6s. Since the 6pi/ 2 
is jcr and | 7 r, there is a greater Pb 7 r population in the ej/ 2 spinors than in the 
nonrelativistic 7 r orbitals. This is offset to some extent by a decrease in the Pb 7r 
population in the e 3 / 2 spinor. There is also a greater 7r population on the metal 
for SnO and GeO. This is reflected in the orbital eigenvalues which are given in 
Table VIII: the 7r eigenvalues are higher by 0.5eV at the DHF level. 

B. RECP results. 

The bond lengths of the monoxides predicted by the RECP calculations are 
nearly all substantially different from the all-electron values. The exceptions are the 
CER+d value for SnO and the CER value for PbO, which are close to the all-electron 
values. These two potentials also gave good results for the hydrides. For PbO the 
CER+d r e value is too long by nearly 0.10A, and the HW value is too short by the 
same amount. Over half the bond lengths have discrepancies with the all-electron 
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values which are greater than 0.03A. The most consistent results were obtained 
with the inclusion of the d shell (in the CER-fd RECP), with the exception for 
PbO just noted. Use of the valence basis from the all-electron calculations in place 
of the supplied basis had little effect on the results except for the SKBJ potentials, 
for which the supplied basis is quite different from the all-electron basis. Even so, 
the basis set changes are not sufficient to bring the results into agreement with the 
all-electron results. 

Deviations from the all-electron dipole moments are also apparent in the 
RECP results. For GeO and SnO, the values are within 0.2 D of the all-electron PT 
value and mostly smaller, but for PbO they are greater by 0.4-0. 6 D. The apparently 
better agreement with the PT value of the HW dipole moment for PbO is due to 
the underestimation of the bond length. Scaling the HW fi e by the ratio of the PT 
and HW r e values yields a value of 6.058, which is in the same range as those of 
the other RECPs. (A calculation with the HW RECP at the PT bond length gives 
a value of 6.386 D). The underestimation of the dipole moment is related to the 
underestimation of the bond length in many cases, and a similar scaling brings the 
values closer to the all-electron result. 6.386 

The harmonic frequencies and IR intensities show similar patterns of de- 
viation from the all-electron results to the stretching mode data of the hydrides. 
The frequencies and intensities are underestimated without the d shell explicitly 
treated in the calculations. The frequency differences are only 10cm —1 for GeO 
but range up to 50cm -1 for PbO. The underestimation of the frequencies in the 
RECP calculations may give a false impression of the importance of the spin-free 
and spin-orbit terms if the results for PbO were compared with experiment and 
with nonrelativistic all-electron calculations, without checking their accuracy by a 
proper calibration of the RECP. At the least, comparisons with PT results should 
be made to determine the range of error incurred in the RECP approximation. 

C. Comparison with other results. 

The RECP results of Igel-Mann et al. [25] on SnO and Basch et al [20] 
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on PbO are generally in good agreement with the all-electron data presented here. 
Although the basis sets are a little different, they are of essentially the same quality. 
Basch et al. use Slater basis functions in their calculations rather than Gaussians. 
The RECP of Igel-Mann et al. is a “full-core” RECP but has a core-polarization 
potential; the RECP of Basch et al. has the 5 d orbital in the valence space, and 
can be labelled a “semi-core” potential. Both potentials thus include the effect of 
the outer core, which is important for a correct description of molecular properties 
[35]. The bond length obtained by Igel-Mann et al. for SnO is a little short and the 
frequency a little high, but the values obtained are better than those obtained from 
the full-core potentials without the core polarization term. The results of Basch et 
al. are in good agreement with the all-electron results, with the exception of the IR 
intensity. This may in part be due to the coarser grid used to fit the dipole function 
around the minimum. 

The minimum Slater basis set bond length of Schwenzer et al. [21] is shorter 
than the present value by 0.04A, as might be expected from a basis without po- 
larization or diffuse functions. The harmonic frequency is in good agreement with 
the present value, but given the discrepancy in the bond length this must be con- 
sidered fortuitous. The ECP results of Datta et al. [22] show close agreement with 
the all-electron values for the bond lengths, both relativistic and nonrelativistic, 
despite the lack of polarization functions, but the harmonic frequencies are much 
too low. Their ECPs do not take the core polarizability into account in any form: 
if they did, the bond length and the frequencies should be larger. The conclusion 
that they draw regarding the basis set deficiency in the calculations of Schwenzer 
et al. [21] is unjustified, since the standard for the evaluation of SCF molecular 
properties should be the basis-set limit all-electron SCF results, not experiment. 


IV. CONCLUSIONS 

Spin-orbit effects are increasingly important as one goes down the periodic 
table, as shown in the calculations on the group IV monoxides. While for GeO, 
PT is adequate to describe the relativistic effects, for SnO spin-orbit effects are 
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sufficiently important to affect the properties noticeably, and for PbO they are so 
important that their omission may cause serious errors in the predicted properties. 
The large spin-orbit interaction in Pb gives rise to bonding between P 1/2 spinors 
on one centre and p 3 / 2 spinors on the other, rather than hybridizing the relativistic 
spinors on each centre to form per and pir bonds with pure spin. 

Several sets of RECPs were calibrated against the all-electron calculations. 
Similar trends to those found for the hydrides were noted in the molecular proper- 
ties, with larger deviations from the all-electron values. None of the sets of RECPs 
gave consistent deviations. The high quality of the CER+d RECP for Sn was con- 
firmed, as was the poor quality of the CER+d and the HW RECPs for Pb. Calibra- 
tion of the results of RECP calculations against equivalent all-electron molecular 
calculations is essential for their use in high accuracy applications, and to avoid 
false conclusions drawn on the basis of comparison with experiment. 
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TABLE I. Bond length r e and relativistic correction to the bond length A rel r e of the XO 
molecules in A. HW: results using Hay and Wadt [5] RECPs. CER: results using RECPs 
from Ref. 7; +d : with outer-core d shell in calculations. SKBJ: results using Stevens et al. 
[6] RECPs. Experimental results are from Huber and Herzberg [19]. 



GeO 

SnO 

PbO 

r e 

HF 

1.597 

1.808 

1.907 

PT 

1.593 

1.797 

1.871 

DHF 

1.594 

1.801 

1.893 

Supplied basis contracted [3s 3p] 

HW 

1.577 

1.739 

1.776 

CER 

1.559 

1.723 

1.893 

CER+d 

1.580 

1.796 

1.967 

SKBJ 

1.560 

1.750 

1.856 

All-electron valence basis 

HW 

1.580 

1.736 

1.776 

CER 

1.561 

1.719 

1.892 

CER+d 

1.581 

1.792 

1.965 

SKBJ 

1.569 

1.746 

1.859 

Other 


1.78 a 

1.864 6 

Expt. 

1.625 

1.833 

1.922 

A rel r e 

Pert. 

-0.0040 

-0.0113 

-0.0363 

DHF 

-0.0030 

-0.0074 

-0.0146 


a Igel-Mann et al., ref. 25. 
b Basch et al., ref. 20. 



TABLE II. Dipole moments fi e and relativistic corrections to the dipole moments A re Ve 
of XO molecules in Debye. HW: results using Hay and Wadt [5] RECPs. CER: results 
using RECPs from Ref. 7; +d: with outer-core d shell in calculations. SKBJ: results using 
Stevens et al. [6] RECPs. Experimental results are from Huber and Herzberg [19]. 



GeO 

SnO 

PbO 

H’e 

HF 

4.076 

5.153 

5.598 

PT 

4.101 

5.160 

5.588 

DHF 

4.107 

5.212 

5.389 

Supplied basis contracted [3s3p] 




HW 

4.079 

5.182 

5.750 

CER 

4.005 

5.003 

6.197 

CER+d 

3.961 

5.124 

6.193 

SKBJ 

4.002 

5.133 

6.044 

All-electron valence basis 




HW 

4.093 

5.073 

5.711 

CER 

4.000 

4.873 

6.160 

CER+d 

3.964 

4.974 

6.151 

SKBJ 

4.082 

5.023 

6.003 

Other 



5.431° 

Expt. 

3.272 

4.32 

4.64 fc 

A rC Ve 




PT 

+0.025 

+0.007 

-0.010 

DHF 

+0.032 

+0.059 

-0.209 


° Evaluated from the SCF dipole moment function of Basch et al., ref. 20. 
b fiQ value. 





TABLE III. Harmonic frequencies in cm -1 and infrared intensities in km mol -1 of XO 
molecules. Intensities are given in parentheses after frequencies. HW: results using Hay and 
Wadt [5] RECPs. CER: results using RECPs from Ref. 7; +d: with outer-core d shell in 
calculations. SKBJ: results using Stevens et al [6] RECPs. Experimental results are from 
Huber and Herzberg [19]. 



GeO 

SnO 

PbO 

NR 

1127(129) 

955(128) 

873(131) 

PT 

1124(134) 

954(136) 

867(174) 

DHF 

1123(133) 

946(136) 

785( 66) 

Supplied basis contracted [3s3p] 

HW 

1113(114) 

928(114) 

819(130) 

CER 

1118(116) 

940(112) 

818(135) 

CER+d 

1137(134) 

949(144) 

834(166) 

SKBJ 

1105(109) 

920(114) 

835(134) 

All-electron valence basis 

HW 

1110(115) 

937(115) 

818(131) 

CER 

1116(120) 

952(112) 

819(135) 

CER+d 

1136(133) 

961(145) 

835(167) 

SKBJ 

1104(121) 

930(115) 

846(137) 

Other 


971“ 

860(105) b 

Expt 

987 

815 

721 


“ Igel-Mann et al., ref. 25. 

b Evaluated from the SCF energy and dipole moment functions of Basch et al., ref. 20. 



TABLE IV. Spinor and orbital eigenvalues in eV of the Ge, Sn, Pb and 0 atoms. 


Ge 


Sn 


Pb 


0 


Relativistic 


Nonrelativistic 


3d 3 / 2 

-44.05 1 
-43.39 . 

3d 

-44.49 

CS 

CO 

4Sj/2 

-15.52 

4s 

-15.16 

4pi/2 

-7.42 i 

4p 

-7.33 

4p 3 /2 

-7.24 . 

4d 3 /2 

-36.36 I 
-35.22 . 

4d 

-37.38 

4^5/2 

5Si/2 

-13.88 

5s 

-13.04 

1/2 

-7.01 i 

5p 

-6.76 

bp 3 /2 

-6.57 . 

bd 3 /2 

-30.99 i 
-28.20 . 

5d 

-33.32 

bds/2 

65j/2 

-15.41 

6s 

-12.48 

6pi/2 

-7.49 i 

6p 

-6.53 

bp 3/2 

—5.99 . 

2sj/2 

-34.08 

2s 

-34.02 

2pi/2 

-16.78 ■ 

2p 

-16.77 

2?3/2 

-16.75 . 



TABLE V. Mulliken population analysis of the valence spinors of XO molecules at r e . The d populations have been omitted, as 
they make little contribution to the valence populations or to the charges. Charges are given relative to the LS'-coupled limit 
(?i/ 2 ) 4 ^ 3 (P 3 / 2 ) 8 ^ 3 for 0, and for Ge, Sn and Pb they are given relative to the jj -coupled limit (pi/ 2 ) 2 - 
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rH 
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CM 

CM 


Pm 


13e 3 / 2 0.14 0.14 0.18 1.82 1.82 1.83 

Charge +0.17 +1.25 -0.46 +0.79 +0.91 +0.01 -0.43 -0.48 -0.90 -0.91 



TABLE VI. DHF integrated spin densities for the valence spinors (Kramers pairs) of the 
XO molecules at r e . The values shown are for the ej/ 2 and e_ 3 / 2 spinors. The time- 
reversed spinors e_j/ 2 and e 3 / 2 have the same densities with opposite spin. The densities 
represent the occupation number of the orbital whose symmetry type given at the head of 
the column. The sum is taken over spin for all occupied valence spinors. 



a{a) 

?r(/3) 

12e ly / 2 

1.000 

0.000 

13 ej / 2 

0.000 

1.000 

14e 1/2 

1.000 

0.000 

5e_ 3 / 2 


1.000 

Sum 

4.000 

4.000 

17*1/2 

0.999 

0.001 

18ej/ 2 

0.015 

0.985 

19 e l/2 

0.984 

0.016 

8 e -3/2 


1.000 

Sum 

3.996 

4.004 

24e!/ 2 

0.992 

0.008 

2hei/2 

0.432 

0.568 

26 e\j2 

0.538 

0.462 


13e_ 3 / 2 

Sum 


1.000 

4.074 


3.926 




TABLE VII. Nonrelativistic Mulliken gross population analysis of the valence orbitals 
of the XO molecules at r e . The d populations have been omitted, as they make little 
contribution to the valence populations or to the charges. 


X 


0 


p Total 


s p Total 


GeO 

9a 

10a 

At: 

Charge 

SnO 

12a 

13a 

6 7T 

Charge 

PbO 

16a 

17a 

97T 

Charge 


0.91 

0.81 

+ 0.18 

0.91 

0.92 

+0.12 

1.10 

0.68 


0.01 

0.57 

0.52 

+ 0.84 

0.02 

0.49 

0.45 

+1.00 

0.01 

0.57 

0.48 

+ 0.89 


0.92 

0.86 

0.63 

+ 0.91 

0.93 

1.42 

0.55 

+ 1.03 

1.12 

1.26 

0.58 

+ 0.96 


0.15 

0.05 

+0.02 

0.08 

0.00 

+ 0.05 

0.06 

0.02 


0.93 

0.56 

3.36 

- 0.90 

0.98 

0.59 

3.44 

- 1.06 

0.82 

0.72 

3.42 

- 1.00 


1.08 

1.12 

3.37 

- 0.91 

1.07 

0.59 

3.46 

- 1.03 

0.88 

0.74 

3.43 

- 0.96 


+ 0.15 


+ 0.05 



TABLE VIII. Spinor and orbital eigenvalues in eV of the XO molecules at r e . 


GeO 


SnO 


PbO 


Relativistic 


12e 1/2 

-16.92 

13e 1/2 

-12.35 -I 

§ e 3/2 

-12.33 . 

14e 1/2 

-11.66 

17e 1/2 

-15.21 

18c 1/2 

-11.01 -I 

& e 3/2 

-11.19 . 

19ei/ 2 

-10.50 

24ej/ 2 

-16.21 

25 e 1/2 

-11.36 -I 

13e 3 / 2 

-10.84 . 

26e!/ 2 

-9.90 


Nonrelativistic 

9a -16.97 

4t r -12.80 

10a -11.69 

12a -14.91 

6tt -11.54 

13a -10.47 

16a -13.96 

97 r —10.61 

17a 


-9.91 
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Abstract 

The ground states of PtH + and PtH 2 and several low-lying states of PtH 
have been studied at the self-consistent field le vel of theory to determine the impor- 
tance of relativistic effects. The results of calculations based on Dirac-Hartree-Fock 
theory, nonrelativistic theory and the spin-free no-pair relativistic approximation of 
Hess are compared to separate the effects of the spin-free terms and the spin-orbit 
terms in the relativistic corrections to the molecular properties. Comparison is 
also made between first-order perturbation theory including the one-electron spin- 
free terms and the Hess method to determine the size of effects beyond first order. 
It is found that the spin-orbit interaction significantly affects the properties and 
energetics of these molecules because of the participation of the Pt 5 d orbital in 
the bonding. Any calculation on Pt compounds which does not explicitly include 
spin-orbit interaction cannot be considered reliable. 


* Mailing address: NASA Ames Research Center, MS RTC 230-3, Moffett Field, 
CA 94035-1000, U.S.A. 


I. INTRODUCTION 


The combustion of hydrogen on platinum surfaces is an important area of 
catalytic chemistry. A number of calculations have been made on models of the dis- 
sociation of H 2 on Pt, treating from one to nine Pt atoms [A1-A16] at various levels 
of theory, most including some treatment of relativistic effects through the use of 
spin-averaged relativistic effective core potentials (RECPs), and most incorporating 
electron correlation effects in some way. It has been recognized that relativistic ef- 
fects must be included in calculations on molecules containing heavy elements, but 
to what lengths it is necessary to go to include them is a question that still needs 
to be answered: for example, in the chemistry of the 6 p block spin-orbit interaction 
is known to be important [A17-A19]. 

The aim of this work is to examine the question of the importance of both 
spin-free and spin-dependent (spin-orbit) relativistic effects for the hydrides of Pt. 
The effect of spin-orbit interaction is assessed by comparison of the DHF results with 
those from the “spin-free no-pair method correct to second-order in the external 
potential” of Hess et al. [A20,A21], hereafter designated PVP. This method is 
derived from the Dirac equation by a transformation decoupling the large and small 
components, and neglecting the spin-orbit terms. Since its Hamiltonian can be used 
in variational calculations, effects beyond first order in a 2 can be obtained. The 
effect of terms beyond first order are obtained by comparison of the PVP results with 
those from first-order perturbation theory (PT) including only the mass- velocity and 
Darwin (MVD) terms [A25,A26]. These are the spin-free relativistic corrections to 
the one-electron operator of order a 2 . The spin- free part of the relativistic correction 
is obtained by comparing the PVP and nonrelativistic results. 


II. THEORY 

The Dirac-Hartree-Fock program described previously [A27,A28] has been 
extended to perform Kramers-restricted open-shell calculations on doublet states, 
where the unpaired electron can be either a single electron outside or a hole inside 


1 


a closed-shell core, or a combination of several such particle or hole states. In the 
presentation of the theory and its implementation, the concepts of Hsu et ah [A29] 
and Faegri and Manne [A30] have been followed. 


For a single electron outside a closed core, the closed- virtual and open- virtual 
Fock operators F cv and F ov are defined by 

p cv pcore _j_ 1_ popen 

jpov m pcore 

with 


( 1 ) 


F core = & + £)<?*, 


k= 1 


(2) 


popen _ q 
1 m ~~ 


where n c is the number of closed-shell Kramers pairs of spinors (the relativistic 
equivalent of orbitals), and 

Gk = Qk + Qk\ Qk = Jk ~ Kk (^) 


with Jk and Kk defined in the usual way but with respect to molecular spinors, and 
the bar indicating time-reversal applied to the specified spinor: spinors \k ) and \k) 
form a Kramers pair. The closed-open Fock operator is obtained from F cv and F ov 
[A30]: 


p co 2 F CV F° v pcore popen 


( 4 ) 


Similarly, for state- averaged doublets, 

pcv pcore _|_ 1 popen 

2 Tl(j ^ 


pov pcore 


pco pcore ^ popen 

where n a is the number of open-shell Kramers pairs, and 

Tic “t” Tlo 


popen ^ ^ popen 

7n=n c - f 1 


( 5 ) 


( 6 ) 
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For a single hole state, the same Fock operators may be used, but it is more conve- 
nient to express them in terms of a hole rather than a particle in the mth Kramers 
pair: 


ptcv pclosed I popen 


jp ov pclosed popen 


p co pclosed 


where p cloaed now includes (? m . For state-averaged doublet holes, 

p cv pclosed 1 popen 


pov _ p 


closed 


p co pclosed 1 popen 


2n 0 -l ' 


( 5 ) 


( 5 ) 


The single Fock matrix method of Faegri and Manne [A30] has been employed in 
the present implementation, with the same choice of diagonal blocks. Hsu et al 
[A29] have argued that in order to ensure that the solution for a core hole state 
corresponds to an excited state, the diagonal core-core and open-open blocks of 
F should be composed of the matrix elements of p clo3ed . This has the effect of 
diagonalizing the Cl matrix of single excitations involving all other core holes, and 
providing some bound to the energy of the core hole state. However, the choice of 
F co already ensures zero interactions with the other core hole states, and the choice 
of P cc and P°° remains arbitrary. 


In the implementation of the open-shell doublet method, only two Fock op- 
erators need to be constructed, F core (or F clo3ed ) and F open . These are constructed 
in the scalar basis, as before, transformed to the 2-spinor basis and combined to 
form a single Fock operator. Advantage is taken in their construction of the fact 
that the same addresses are required in the gather /scatter operation for both Fock 
operators. In the present calculations on PtH, the extra time taken for the second 
Fock operator amounted to less than 20% of the time for the first Fock operator. 
The scheme for calculating the addresses has also been improved and time-reversal 
symmetry implemented in the construction of F, resulting in a saving in CPU time 
of more than 30%. 
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III. COMPUTATIONAL DETAILS 


Moderately large Gaussian basis sets have been used for Pt and H, so that 
the results should be close to the DHF limit. For Pt, the 22sl6pl3d8/ primitive 
basis set of Faegri [A31] was taken as a starting set. An extra s function was added 
to the Valence space and the three outer s functions reoptimized. Three p and 
two d functions were added in an even-tempered series with ratio 0.4, and two / 
functions with exponents 0.2 and 0.08 were added for polarization of the 5 d orbital. 
In the relativistic basis set an extra tight p function with exponent 6.8 x 10 s was 
added to make up part of the deficiency in the p-^j 2 space [A32]. The addition 
of this function lowered the DHF total energy by 0.5 Eh, leaving a discrepancy 
of 0.152 Eh between it and the DHF limiting value, compared with 0.01 Eh for 
the nonrelativistic discrepancy. The core was generally contracted to 5s4p3dl/ 
using the atomic SCF orbitals, leaving a primitive valence set of 4s4p4d2/. For 
hydrogen a 3s 2 p contracted set was used, taking the 5s primitive set of Huzinaga 
[A33] contracted to 3s, and two p functions with exponents 1.25 and 0.45 [A34]. 
The relativistic contraction coefficients were obtained from an adaption of GRASP 
[A35] described previously [A28]. 

The results from nonrelativistic (NR), first-order perturbation theory (PT) 
and the PVP method were obtained with the MOLECULE/SWEDEN [A36] pack- 
age. DHF results were obtained with the program described above. The value of 
the speed of light in atomic units was taken to be 137.0359895. SCF calculations 
for PtH and PtH + were carried out at 5 internuclear separations around the mini- 
mum distance with a spacing of O.lA. The energies were fit to a quartic function 
in the internuclear separation. SCF calculations for PtH 2 were carried out at 4 
Pt-H distances around the minimum distance with a spacing of 0.07A and at 5 
H— Pt-H angles around the minimum spaced at 2°. The energies were fit to a cubic 
function in the two coordinates. Infrared (ir) data were obtained with the program 
INTDER [A37] All calculations were performed on the Computational Chemistry 
CONVEX C-210 and Central Computing Facility CRAY Y-MP/864 computers at 
NASA Ames Research Center. 
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IV. RESULTS AND DISCUSSION 


A. Pt atom 

It has been pointed out by a number of authors [A16, A38] that in order to 
obtain the correct ordering of molecular states, it is necessary to have the correct 
ordering of the atomic states. This is particularly so for transition metal atoms 
where there are three configurations that contribute low-lying states to the atomic 
spectrum, viz. d n , d n ~ 1 s and d n ~ 2 s 2 . These remarks have been made in the context 
of the inclusion of correlation effects. As well as correlation effects, for heavy atoms 
such as Pt, relativistic effects are of such great importance that a correct description 
of the atomic states cannot omit them. Many calculations which include relativistic 
effects via RECPs do not include spin-orbit effects, since these require considerable 
extra work. The comparisons with experiment made in these studies average over 
the J components of an LS term. A cursory examination of the levels of the d 9 s 
configuration from the optical spectrum of Pt [A39] shows that electrostatic and 
spin-orbit interactions are of comparable magnitude. A correct description of the 
Pt atom will therefore have to take spin-orbit effects into account, and there is 
some doubt whether averaging over J values will give a meaningful term energy 
with which to compare, since the levels of a given LS multiplet do not obey the 
Lande interval rule. 

The results of DHF calculations using GRASP [A35] on the (d+ s) 10 config- 
uration complex, followed by Cl within the complex, are presented in Table I. The 
state designations axe given in both jj and LS coupling. The effect of the Breit 
interaction and estimated quantum electrodynamic (QED) effects are also included. 
The strong configuration mixing makes the assignment of states to given terms dif- 
ficult, and the multiplet designations given in the experimental spectrum may have 
to be revised. The J quantum number is the only reliable designator, and based 
on this, the calculated levels agree fairly well with the experimental levels, when 
the designation of the d 8 s 2 3 P 2 and the d 9 s 1 1 D 2 have been reversed. These two 
states are stongly mixed in any case, so that the labelling is somewhat academic. 
The results are sensitive to the weights of the three configurations used in the SCF 
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optimization, with the best overall results (given here) obtained with equal weights 
of the d s s 2 and d 9 s configurations, and 10% relative weight of the d 10 configura- 
tion. However, use of the orbitals from an SCF calculation on the d 9 s configuration 
lowered the energy of the d 10 state by 10000 cm -1 while raising the d s s 2 state ener- 
gies by 5000 cm -1 ; using 1000% relative weight of the d 10 configuration gave a d 10 
level energy of 5889 cm -1 , which is close to the experimented value of 6140 cm -1 . 
Clearly, a single orbital set cannot adequately describe all three configurations, and 
correlation effects, including orbital relaxation, must be included to obtain a more 
quantitative prediction of the energy levels. It is noteworthy that the inclusion 
of the Breit interaction made relative shifts of energy levels of up to 800 cm -1 , 
indicating that it must also be included for quantitative (“chemical”) accuracy. 

B. PtH+ 

Little is known about PtH + experimentally. The ground state is asserted to 
be 3 A by Wang and Pitzer [A3], which can have w values of 1, 2 and 3. Balasub- 
r am an i art ’s [A9] calculations at 3.1a<), which include both spin-orbit and electron 
correlation effects, show that the ground state has 0 + symmetry and is dominated 
by the nonrelativistic 3 E + configuration, but with a substantial contribution from 
the 3 H; the 3 A 3 is the next state, more than 4000 cm -1 higher in energy. 

DHF calculations have been performed on the lowest 0 + state of PtH + , and 
NR, PVP and PT calculations on the 1 E + , 3 n and 3 A states. The results, including 
basis set superposition error (BSSE) calculated by the counterpoise method [A40], 
axe reported in Table II. The fragments for the interaction energies were taken to 
be Pt + and H, justified by the ionization potentials (IPs) of H and Pt and by the 
Mulliken populations, given in Table III at r = 3.1ao, which place most of the 
molecular charge on the Pt. Only the PT calculations were significantly affected by 
BSSE. The largest effect was a 20cm -1 increase in the 1 E + harmonic frequency, and 
the effects on D e were less than 0.2kcal mol -1 . At the NR ROHF level, the 3 H and 
the 3 A curves were repulsive in the vicinity of the 1 E + minimum, but the 3 A curve 
was less repulsive than the 3 H curve. Addition of the mass- velocity and Darwin 
(MVD) terms stabilized all three states, but the 3 H state still only had a shallow 
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minimum, at a value of r greater than 3.5ao< The stabilization can be related to the 
s population on Pt, from which the principal contribution of the MVD operator to 
the valence energy comes, and which is largest for the 3 A state and smallest for the 
1 E + state. The stability of the states correlates with the amount of charge transfer 
from H to the Pt + ion. Electron correlation effects undoubtedly lower the 1 E + state 
more than either of the triplet states, and from the current calculations and those 
of Balasubramanian it can be concluded that the ground state of PtH + is indeed 
the 1 E + (0+) state. 

The changes to the bonding orbitals induced by relativity are illustrated 
by the DHF Mulliken gross populations at r = 3.1ao, also given in Table III. 
The expansion of the Pt 5 d orbital enables charge to be transferred from it to the 
relativistically stabilized Pt 6s and the H Is orbitals. The latter transfer reduces 
the charge on H, making the molecule more like Pt + — H. The strong influence of 
spin-orbit interaction is shown by the breakdown of the DHF populations by spin- 
orbit components and spin densities, displayed in Table IV. Both valence spinors 
have almost equal contributions from a and f3 spin, a situation about as fax from A-s 
coupling as possible. The 20e!/ 2 spinor has almost no Pt 5 d 5 / 2 contribution, and 
the contributions to the 21ej/ 2 spinor from the two 5 d spin-orbit components are by 
no means in the X-s coupling ratio. The dominance in the MSs of one 5 d spin-orbit 
component on Pt has an inductive effect on the H p spinors, forcing one of them 
to give a negligible contribution. This effect has been noted previously [A19], and 
occurs because the combination of spinors with opposite k* on two centres gives 
bonding functions in both a and (3 spin spaces, whereas combinations with the same 
k give a bonding function for one spin with an antibonding functions for the other 
spin. Table IV also shows that the charge on the Pt atom is due to removal of an 
electron from a 5 d 5 / 2 spinor, with the rest of the charge transfer from the 5 d shell 
distributed between both spin-orbit components. 

The effect of spin-orbit interaction on the molecular properties can readily 

* Quantum number of the relativistic angular operator K = —1 — 2S • L, which is 
related to the spin-orbit interaction. 
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be explained by the mixing of the 3 Hq+ state into the 1 S^ + state. The 3 IIo+ state 
has a larger bond length and smaller frequency, and the interaction between the 
two will increase D e . The importance of spin-orbit interaction is shown by the 
size of the differences between the PT results and the DHF results: 0.02A in r e , 
120cm -1 in u e and 7.5kcal mol -1 in D e . The extent of the mixing can be judged 
from the integrated spin densities, which show a a — > 7 r promotion of 0.23e. Clearly, 
spin-orbit effects cannot be ignored in this molecule. 

It is necessary to exercise some degree of caution with the results of PT 
calculations involving the MVD terms. It has been shown for some transition metal 
compounds that the results of PT calculations are not stable with respect to the 
contraction scheme [A41] . A systematic study of the effect of core orbital relaxation 
on the d 10 l S -+ d 9 s 3 D excitation energy of Pd was made by Blomberg and 
Wahlgren [A42], in which they showed that the relaxation of the outer core orbitals 
was necessary to obtain consistent PT results. In contrast, neither the NR SCF 
nor the no-pair approximation of Almlof et al. [A43] showed much dependence of 
the excitation energy on the core relaxation. A series of calculations was performed 
to assess the effect of contraction on the properties of PtH + at the NR and PT 
levels of theory. The valence basis was increased by uncontracting the outermost 
primitive functions from the core, one at a time, until the properties stabilized. 
When the second s and p functions were uncontracted the core parts of the 5s and 
bp were discarded for reasons of linear dependence. The results are given in Table 
V. The most dramatic change in the PT properties occurred with the uncontraction 
of the first s f un ction — a not unreasonable result, considering the nature of the 
mass-velocity and Darwin operators. Given the changes, the PT results reported in 
Table II are for the bsbp5d valence basis. The very small changes in the NR values, 
coupled with the experience of Blomberg and Wahlgren, give grounds for confidence 
that the DHF results are stable. Gropen et al. [A16], in their calculations on PtH, 
have also found that a flexible valence basis is necessary when using PT. 

C. PtH 

Spin-orbit effects are dominant for this molecule, and the structure of the 
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low-lying states cannot be understood without it. In A -s coupling there are three 
terms that arise from the d hole on Pt, 2 A, 2 13 and 2 E + . The first two are split 
by spin-orbit effects to give 2 A 5 / 2 > 2 A 3 / 2 > 2 Il 3 / 2 ? and 2 Hi /2 states. The 2 A 3 /2 and 
2 n 3/2 states are mixed by the spin-orbit interaction as are the 2 Hi/ 2 and 2 E]y 2 . In 
u}-u) coupling there are two sets of states, arising from the 5<f 5 /2 and 5 ^ 3/2 holes on 
Pt, the former having u> = 5/2, 3/2 and 1/2 and the latter having u> = 3/2 and 
1/2, the former correlating with the atomic states J = 3, 2*, designated 

in LSJ coupling as 3 D^ and 3 Di, and the latter correlating with the atomic states 
J = 2, 1, designated in LS coupling as l Di and 3 D\. The results of three 
sets of calculations are presented in Table VI, a 5-state average DHF calculation, a 
3-state average DHF calculation and separate calculations on the the lowest three 
states. The 5-state average calculation clearly shows the grouping of the molecular 
states based on the jj-coupled hole states: the two groups are separated by the 5 d 
spin-orbit splitting of about 10000 cm -1 . The importance of the spin-orbit splitting 
is also reflected in the molecular properties. The 5/2 and 3/2(1) states, which in A-s 
coupling would be spin-orbit components of the 2 A state, have significantly different 
bond lengths, harmonic frequencies and dipole moments. The differences are too 
large for a small perturbation. There is more similarity between the 3/2(1) and 
1/2(1) states, which are not components of a spin-orbit doublet. The differences 
in the molecular properties also cast doubt on any procedure which includes the 
spin-orbit splitting as a constant shift applied to the non-relativistic curves (or even 
curves obtained with RECPs) on the basis that the splitting is local to the heavy 
atom. The involvement of the spin-orbit split orbital in bonding undermines the 
assumption of locality. 

The results of the NR and PT calculations are presented in Table VII. Similar 
problems with the PT calculations were encountered as for PtH + . Results for both 
the 4s4p4d and 5s5p5d valence basis sets are reported. The seriousness of the errors 
incurred with the more highly contracted basis is emphasized by the fact that the 
predicted ground state changes when the basis is uncontracted. The importance 
of relativistic effects in stabilizing the molecular states is obvious. At the NR SCF 

* The square brackets indicate a hole in a filled shell 
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level, only the 2 E + state is bound. However, comparison of the PT results with 
the DHF results shows that the omission of the spin-orbit interaction is at least as 
serious as the omission of the non-fine-structure terms. The spread of states is only 
3000cm -1 , which is similar to the spread of states based on the 5d 5 /2 hole in the 
DHF calculations. Thus the valence electrostatic interactions are only a third of 
the size of the spin-orbit splitting. 

The Mulliken gross populations are presented in Tables VIII and IX. The 
most significant difference between the NR and DHF populations is the transfer of 
charge from the d shell into the s shell. It is noteworthy that the atomic charges of 
the 5/2 and 3/2 states is lower than that of the 2 A and 2 H states, but the 1/2 state 
has the same charges as the 2 E + state. 

Further evidence for the importance of spin-orbit effects in determining the 
state designations is obtained from the atomic spinor contributions to the Mulhken 
populations and the integrated spin densities of the valence MSs. These are pre- 
sented in Tables X and XI. If the molecule is well-described in the X-s coupling 
scheme, the ratio of the j = i + 1/2 and j = t — 1/2 component densities should 
be ^ + m + l : i — m for spinors with u> = m + 1/2 for pure a spin, and the reverse 
for pure /? spin. Few of the spinors have ratios that are close to these values. The 
integrated spin densities are perhaps a better indication of the mixing. The 20ej/2 
and 22e 1 ^ 2 spinors are mostly a and the 21 ej /2 mostly 7 r, but with strong mixing, 
especially for the 1/2 state. The sums of the densities give the orbital occupation 
numbers. While the 5/2 state is very close to a pure 5 hole state, the 3/2 state 
has a strong mixing of the 7 r and delta holes, and the 1/2 state a smaller but still 
significant mixing of the <x and 7 r holes. 

There has been a number of previous calculations on PtH [A1,A3,A8,A13,A16], 
mainly in connection with the 2 A — 2 separation. Wang and Pitzer [A3] report 
RECP SCF calculations, and examine the effect of spin-orbit interaction via Cl and 
correlation via singles and doubles (SD)CI calculations. With the spin-orbit inter- 
action but without correlation, the 1/2(1) and 3/2(1) states have excitation energies 
of -1229 cm -1 and 4082 cm -1 relative to the 5/2 state; correlation changes these 
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to 1008 cm -1 and 2742 cm -1 . The 1/2(2) state is 11810 cm -1 above the 1/2(1) 
state, but the 3/2(2) state is 19350 cm -1 higher than the 3/2(1) state. These 
calculations, while they demonstrate the importance of the spin-orbit interaction, 
clearly have some deficiencies. Rohlfing et al. have performed a systematic set of 
RECP calculations on the states of NiH, PdH and PtH, incorporating correlation by 
M0ller-Plesset (MP) perturbation theory up to fourth order. Their results demon- 
strate the importance of inclusion of the semi-core 5s and 5 p orbitals on Pt in the 
valence space, particularly when electron correlation is included. The correlated 
results show a remarkably small change in the excitation energies of the 2 A and 
2 n states. Their UHF results have been included in Table ? for comparison with 
the present ROHF results. The properties display significant deviations from the 
PT and the PVP results — shorter bond lengths, higher harmonic frequencies and 
dissociation energies. Although their valence basis set is somewhat smaller, it does 
contain an / function on Pt and a p function on H, and should not give results too 
different from the basis used in this work. If the contraction of the bond length 
due to correlation for the 2 A state obtained by them is applied to the DHF results, 
the bond length 5/2 state is in very good agreement with the experimental value of 
1.5285A. Nevertheless, the results of Rohlfing et al. [A8] for T e using a simple model 
of the spin-orbit splitting are in reasonable agreement with the present values. 


D. PtH 2 

The equilibrium geometries, dipole moments and energies for dissociation 
into Pt 1 So and H 2 are presented in Table XII. No nonrelativistic (NR) SCF results 
are given because at this level, PtH 2 is unbound — except perhaps for a van der 
Waals minimum between the Pt atom and the H 2 molecule. This is not surprising, 
since the nonrelativistic ground configuration of Pt atom is d 10 , a closed-shell config- 
uration. The energy required to break up a closed-shell configuration by promotion 
into the next unoccupied orbital is usually too high for the formation of strong 
bonds, and such atoms tend to form van der Waals complexes. For Pt, the d — s 
promotion energy is considerably reduced by the lowest order relativistic effects, so 
that the inclusion of PT makes PtH 2 stable with respect to both the d 10 and the 
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ground state d 9 s 1 S D Pt atom asymptotes, in the latter case by 3.5kcal mol 1 . 

What is required is promotion of an electron from the d orbital into the 
empty s orbital [A5,A6]. Relativistic effects stabilize the s orbital and destabilize 
the d orbital sufficiently for Pt that the ground configuration is cPs 1 . The s or- 
bital stabilization is already obtained with TT, and the molecule is stabilized by 
15kcal mol -1 relative to the d 10 asymptote, and is still bound by 3.5kcal mol -1 
relative to the ground d 9 s 3 3 D state of Pt and H 2 . 

There is some little confusion over the issue of the bonding in PtH^. In 
his conclusions, Balasubramanian [A9] predicts that “the Pt( 1 5'o) state would in- 
sert into H 2 almost spontaneously to form the bent PtH 2 ground state while the 
3 Z >3 state would not react with H 2 spontaneously”. These conclusions are based on 
RECP/CASSCF/CI calculations which show only a small barrier in the 1 A\ surface 
to the insertion, but a large barrier in the 3 Ai surface, with a linear equilibrium ge- 
ometry and a higher energy than the asymptote. The Mulliken populations he gives 
indicate a substantial Pt 6s population (0.75e), which suggests that d — s promo- 
tion is necessary for bonding, a conclusion arrived at by Low and Goddard [A5,A6]. 
Since the ground state of PtH 2 is (in LS coupling) a singlet, its formation must 
occur from a singlet state of the atom. When spin-orbit interaction is introduced, 
this requirement is removed. However, dissociation of PtHj at the Dirac-Fock level 
must still occur to a state of the atom in which all electrons are paired, and the 
lowest such state is the d 10 state. As the dissociation limit is approached, the curve 
for this state must undergo avoided crossings with triplet states arising from the 
d 9 s 3 3 Dz, d 9 s l 3 £> 2 and d 8 s 2 3 F^ asymptotes, all of which can contribute a molecular 
state with Ay double-group symmetry. 

Spin-orbit effects have a marked impact on the geometry. The DHF bond 
angle is nearly 5° greater than the PT bond angle. This is much larger than the 
difference for PbH 2 [A18], which was only 1°. The analysis of the electron density 
into contributions from the single-group irreps gives insight into this increase in bond 
angle. In the valence region, electron density is transferred by spin-orbit effects into 
the by component from all the others, but mainly from the ay component. The 
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greater density in the 61 component increases the repulsion between the bonding 
pairs which increases the bond angle. 

Such an increase in bond angle due to spin-orbit effects has important im- 
plications for the dissociation of H 2 on Pt surfaces. Sj0voll’s [A15] calculations on 
Pt 9 H 2 indicate that the dissociated H 2 molecule in fourfold hollow sites on the 100 
surface is less stable than the locally PtH 2 -like configuration. Spin-orbit effects, 
in opening up the bond angle, would then facilitate the dissociation of H 2 into a 
chemisorbed state. Moreover, the lowering of symmetry due to the spin-orbit inter- 
action may remove the necessity for preparing the cluster for bonding to H 2 : in (7 2 v 
symmetry there is only one fermion double-group irreducible representation (irrep), 
so that the MSs can freely change their composition in terms of single-group basis 
functions. 

The bond angle obtained in the PT calculations is similar to the value 
obtained by Sj0voll [A15], but smaller than that obtained by Balasubramanian 
[A9]. Both these calculations employed RECPs and included correlation via 
CASSCF/MRCI calculations. Sj0voll’s basis included / functions and had the 4 d, 
5s and 5 p semi-core orbitals frozen, where Balasubramanian’s had no / functions 
and no frozen semi-core orbitals. 


IV. CONCLUSIONS 

Strong spin-orbit effects are evident in the hydrides of Pt. The grouping 
of the low-lying states of PtH is determined by the large Pt atomic 5 d spin-orbit 
splitting, and the potential energy curves are not parallel, as might be expected if 
spin-orbit effects were small. The geometries and harmonic frequencies of the states 
PtH"*" and PtH are significantly altered by spin-orbit effects; the bond tingle in PtH 2 
is increased by nearly 5° by spin-orbit effects. While correlation effects are very 
important for transition metal compounds, the spin-orbit splitting is of comparable 
importance. Whether it can be treated as a perturbation in any accurate calculation 
on Pt compounds is highly questionable. 
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TABLE II. Properties of the low-lying states of PtH + from NR, PT and DHF calculations 
on each individual state. Bond lengths r e are given in A, harmonic frequencies w e and 
excitation energies T e in cm -1 , dissociation energies D e in kcal mol . 


State 

T e 

CJe 

T e 

D e 

NR 

2 E+ 

1.497 

2464 

0 

9.7 

PT 

1 E + 

1.501 

2546 

0 

20.6 

3 A 

DHF 

1.516 

2550 

1160 

17.3 

0+ 

1.522 

2426 

0 

28.1 


TABLE III. Gross Mulliken populations and charges for the NR 1 E + , 3 II and 3 A states and 
the DHF 0 + state of PtH + at r = 3.1a 0 . The populations of the valence a orbitals and e 1 / 2 
spinors are displayed for each state, as well as the total populations and the net charge. 
The / populations on Pt have been omitted, as they make little contribution to the valence 
populations or to the charges. The charges for Pt are given relative to the d 10 state of the 
neutral atom. 




Pt s 

Pt p 

Pt d 

Total 

H s 

H p 

Total 

3 E+ 

13cr 

0.074 

0.039 

1.069 

1.186 

0.794 

0.020 

0.814 


Charge 

-0.074 

-0.032 

+0.944 

+0.834 

+0.200 

-0.034 

+0.166 

3 n 


13a 

0.003 

0.002 

1.715 

1.720 

0.262 

0.018 

0.280 


140- 

0.130 

0.069 

0.108 

0.312 

0.689 

-0.002 

0.687 


Charge 

-0.131 

-0.058 

+1.187 

+0.993 

+0.037 

-0.030 

+0.007 

3 A 


13er 

0.009 

0.015 

1.654 

1.680 

0.307 

0.013 

0.320 


140- 

0.178 

0.057 

0.139 

0.376 

0.625 

-0.002 

0.624 


Charge 

-0.187 

-0.067 

+1.215 

+0.958 

+0.063 

-0.021 

+0.042 

0+ 


20e 1 / 2 

0.036 

0.015 

1.671 

1.723 

0.265 

0.012 

0.277 


21e 1/2 

0.292 

0.054 

0.992 

1.341 

0.646 

0.013 

0.659 


Charge 

-0.326 

-0.066 

+1.344 

+0.948 

+0.084 

-0.032 

+0.052 



TABLE IV. DHF gross Mulliken populations and integrated spin densities for the spin-orbit 
components of the valence spinors, total charges and total integrated spin densities of the 
0 + state of PtH + at r = 3.1a 0 . The f populations on Pt have been omitted, as they make 
little contribution to the valence populations or to the charges. The charges for Pt are given 
relative to the d 10 state of the neutral atom. 


Pt H Spin 



P 1/2 

P 3/2 

^ 3/2 

d 5 / 2 

P 1/2 

P 3/2 

a 

P 

20ei/ 2 

0.004 

0.011 

1.667 

0.004 

0.000 

0.013 

0.540 

0.460 

21e!/ 2 

0.026 

0.028 

0.219 

0.772 

0.012 

0.002 

0.576 

0.424 

Charge 

-0.030 

-0.036 

+0.115 

+1.230 

-0.011 

-0.021 







TABLE VI. Properties of the low-lying states of PtH from DHF calculations. The states are 
designated by D with the energetic order index in parentheses. Three sets of calculations 
are described: a 5-state average, a 3-state average and a set of separate calculations on each 
state. Bond lengths r e are given in A, harmonic frequencies u; e and excitation energies T e 
in cm -1 , dipole moments in Debye with polarity Pt + H _ . 


State 

r e 

CJ e 

T e 

fl e 

5-state average 

5/2(1) 

1.551 

2241 

0 


3/2(1) 

1.587 

2074 

3226 


1/2(1) 

1.590 

2044 

3603 


3/2(2) 

1.575 

2166 

11290 


1/2(2) 

1.588 

2129 

11904 


3-state average 

5/2(1) 

1.552 

2238 

0 


3/2(1) 

1.587 

2078 

3226 


1/2(1) 

1.583 

2085 

2882 


Separate 

5/2(1) 

1.551 

2234 

0 

2.316 

3/2(1) 

1.584 

2080 

3365 

2.558 

1/2(1) 

1.573 

2094 

2712 

2.463 


TABLE VII. Properties of the low-lying states of PtH from NR and PT calculations on 
each individual state. Two PT calculations are shown, with different valence basis sets. 
Bond lengths r e are given in A, harmonic frequencies u e and excitation energies T e in cm -1 , 
dissociation energies D e in kcal mol _1 and dipole moments p e in Debye. The NR D e is given 
relative to Pt d 10 ; the PT D e is given relative to Pt d 9 s. RHM: ECP(nc/18) UHF results 
of Rohlfing et al., ref A8. 


State 

r e 

u; e 

Te 

D e 

jLX e 

NR 

2 S + 

1.652 

1754 

0 

14.2 

2.919 

PT 4s4p4d 
2 A 

1.579 

2035 

0 

38.8 


2 n 

1.655 

1815 

3229 

29.6 


2 S + 

1.575 

1960 

650 

36.9 


PT 5s5p5d 

2 A 

1.593 

2014 

246 

37.9 

2.932 

2 n 

1.670 

1811 

3333 

29.1 

3.472 

2 S + 

1.576 

1949 

0 

38.6 

3.234 

RHM 

2 A 

1.539 

2292 

0 

48.2 


2 n 

1.597 

2050 

4597 

35.0 


2 S + 

1.551 

2194 

1371 

44.5 



TABLE VIII. Mulliken population analysis from DHF calculations on the lowest w = 5/2, 
3/2 and 1/2 states of PtH at r = 3.0a o . The populations of the valence e^/i spinors 
are displayed for each state, as well as the total populations and the net charge. The 
/ populations on Pt have been omitted, as they make little contribution to the valence 
populations or to the charges. The atomic charges for Pt are given relative to the d 9 s state. 


Pt 


H 


5/2 


3/2 


1/2 


s p d Total s p Total 


20ei/ 2 

0.045 

21ei/ 2 

0.033 

22e!/ 2 

1.024 

Charge 

-0.100 

20e ly /2 

0.042 

21^1/2 

0.037 

22e!/ 2 

0.969 

Charge 

-0.045 

20e 1 / 2 

0.107 

21^1/2 

0.305 

22e 1/2 

0.309 

Charge 

+0.282 


0.034 

1.332 

0.019 

1.727 

0.004 

0.594 

-0.053 

+0.355 

0.024 

1.347 

0.018 

1.701 

0.014 

0.624 

-0.047 

+0.334 

0.024 

1.236 

0.034 

1.129 

0.000 

0.644 

-0.053 

+0.001 


1.411 

0.576 

1.780 

0.212 

1.625 

0.375 

+0.199 

-0.171 

1.411 

0.575 

1.756 

0.234 

1.612 

0.388 

+0.240 

-0.208 

1.367 

0.620 

1.471 

0.522 

0.952 

0.047 

+0.227 

-0.197 


0.012 

0.589 

0.008 

0.220 

0.000 

0.375 

-0.028 

-0.199 

0.014 

0.589 

0.010 

0.244 

0.000 

0.388 

-0.032 

-0.240 

0.013 

0.633 

0.007 

0.529 

0.001 

0.048 

-0.031 

-0.227 



TABLE X. Mulliken population analysis of the atomic spin-orbit components of the valence 
spinors for the lowest u> = 5/2, 3/2 and 1/2 states of PtH at r = 3.0a o . The / populations 
on Pt have been omitted, as they make little contribution to the valence populations or to 
the charges. The charges for Pt are given relative to the d^ 2 d^ 2 s^ 2 state of the atom. 


Pt H 




P 1/2 

P 3/2 

d 3 / 2 

^ 5/2 

P 1/2 

P 3/2 

5/2 


20e 1/2 

0.013 

0.021 

1.142 

0.190 

0.001 

0.011 


21ei/ 2 

0.013 

0.006 

0.656 

1.071 

0.009 

-0.001 


22 e !/ 2 

0.003 

0.001 

0.104 

0.491 

0.000 

0.000 


Charge 

-0.028 

-0.025 

+0.100 

+0.255 

-0.009 

-0.019 

3/2 


20e 1 / 2 

0.009 

0.015 

1.120 

0.227 

0.001 

0.013 


21 ^ 1/2 

0.014 

0.004 

0.718 

0.983 

0.011 

-0.001 



0.007 

0.008 

0.078 

0.546 

0.000 

0.000 


Charge 

-0.028 

-0.019 

+0.141 

+0.193 

-0.013 

-0.019 

1/2 


20 e !/ 2 

0.008 

0.017 

1.200 

0.036 

0.000 

0.013 


21 ej /2 

0.022 

0.012 

0.677 

0.452 

0.008 

-0.001 


22e 1/2 

-0.001 

0.000 

0.032 

1.256 

0.002 

0.000 


Charge 

-0.028 

-0.025 

+0.109 

-0.108 

-0.009 

-0.021 


TABLE XI. DHF integrated spin densities for the valence spinors (Kramers pairs) of the 
lowest ijj = 5/2, 3/2 and 1/2 states of Pt H at r — 3.0a o . The values shown are for the 
e i/2) e -3/2 and e 5/2 spinors. The time-reversed spinors e_!/ 2 , £ 3/2 an d e_ 5 / 2 have the same 
densities with opposite spin. The densities represent the occupation number of the orbital 
whose symmetry type given at the head of the column. The sum is taken over spin for all 
occupied valence spinors. (The core sums give the expected integer occupation numbers.) 


5/2 


3/2 


1/2 


<r{a) tt(/?) 6(a) 


20e 1/2 

0.847 

21ei/ 2 

0.175 

22ej/ 2 

0.981 

H e -3/2 

12e_ 3 / 2 

5e5/ 2 

Sum 

4.006 

20ej/ 2 

0.865 

21^1/2 

0.172 

22ej/ 2 

0.968 

H e -3/2 

12e_ 3 / 2 

5e5/ 2 

Sum 

4.010 

20ej/ 2 

0.739 

21ei/ 2 

0.430 

22ej/ 2 

0.836 

H e -3/2 


12e_3/ 2 

5^5/2 

Sum 

3.174 


0.153 


0.825 


0.019 


0.875 

0.125 

0.125 

0.875 


1.000 

3.994 

3.000 

0.135 


0.828 


0.032 


0.584 

0.416 

0.416 

0.584 


1.000 

3.574 

3.416 

0.261 


0.570 


0.164 


0.553 

0.447 

0.447 

0.553 


1.000 

3.826 

4.000 


TABLE XII. Bond length, r e , in A, bond angle, 0 e , in degrees, dipole moment, /i e , in Debye 
and energy of decomposition into Pt and H 2 , A E, in kcal/mol into Pt and H 2 of PtH 2 - 


0 e He A e 


PT 

1.534 

83.3 


13.3 

PVP 

DHF 

1.525 

88.1 

2.503 

35.1 



fir !<=^\ 


SUMMARY 


The investigation into the appearance of intmder states from the neg- 
ative continuum then some of the t.o-electron integrals involving the small 
component of the uave function were omitted has been completed and published 
[1] . The work shows that provided all integrals involving cor. contracted 
functions in an atomic general contraction are included, or that the core 
functions are radially localized, meaningful results are obtained and in- 
truder states do not appear. Reprints of the paper are attached. 


In the area of program development , the Dirac-Hartree-Fock (DHF) pro- 
gram for closed-shell polyatomic molecules has been extended to permit 
Kramer s-restricted open-shell DHF calculations with one electron in an open 


shell or one hole in a closed shell, or state-averaged DHF calculations over 
several particle or hole doublet states. Code for two open shells is still 
being tested, as is code for limited MCSCF, for situations where the ground 
state at the DHF level cannot be described by a single determinant, but is a 


linear combination of a few determinants related by double excitation, from 
one Kramers pair into another. One application of the open-shell code was to 
the KO molecule [2] . Preprints of this paper are attached. 


Another major area of program development which is under development is 
the transformation of integrals from the scalar basis in which they are gener- 
ated to the 2-spinor basis employed in parts of the DHF program, and thence to 
supermatrix form. The reason for these development, is that the DHF program, 
while written efficiently for the circumstances under which it was expected 
to be used, is not the most efficient implementation possible. The code was 




originally developed under the assumption that disk space would not be avail 
able to sort and transform the integrals, and therefore the unordered inte- 
grals would have to be used in the construction of che Fock matrix. With the 
experience gained over the past year, particularly concerning the omission 
of small component integrals, and with the increase in availability of disk 
space, it is now possible to consider transforming the integrals. The use of 
ordered integrals, either in the scalar basis or in the 2-spinor basis, would 
considerably speed up the construction of the Fock matrix, and even more so if 
supermatrices were constructed. Furthermore, in order to proceed beyond the 
SCF level and include electron correlation it is necessary to transform the 
integrals to the molecular 4-spinor basis. Therefore, a considerable amount 
of effort has been spent on analyzing the integral ordering and transformation 
for the DHF problem. Much of this work was used in preparation for the NATO Ad- 
vanced Summer Institute in Vancouver, BC, in August 1992. 


The work of assessing the reliability of the relativistic effective core 
potentials (RECPs) available in the literature has been continued with calcu- 
lations on the group IV monoxides . The perturbation of the metal atom provided 
by oxygen is expected to be larger than that provided by hydrogen and thus pro- 
vide a better test of the quality of the RECPs. The results of this study have 
been submitted for publication, and preprints are attached [3] . 

Calculations on the platinum hydrides PtH, PtH + and PtH 2 have been car- 
ried out at the nonrelativistic (NR), perturbation theory (PT) and DHF levels. 
The DHF calculations employed the new open-shell code described above. Geome- 
tries, dipole moments, harmonic frequencies, infrared intensities and disso- 




ciation energies have been calculated. For purposes of separating spin-orbit 
effects from the non-fine-structure effects and for determining the validity 
of the first-order PT approximation, it is proposed to include in the results 
for the platinum hydrides calculations with the spin-free no-pair code of Hess 
[4] . A preliminary draft of the paper without these calculations is attached. 
These results were presented at the West Coast Theoretical Chemistry Confer- 
ence, held in May 1992 at Pacific Northwest Laboratory. 
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